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DEFINITION OF AXIS SYSTEMS 


Unless otherwise Indicated, all moment coefficients and derivatives 
are referenced to the fuselage-centerline-oriented body system of axes. 
All longitudinal force coefficients and derivatives are referenced to 
the flight path (or wind) system of axes. Thene axis systems are de¬ 
fined in the sketch below. 
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SECTION I 


introduction 

A. BACKGROUND 

For modern military aircraft, manageable flight at high, near-stall 
angles of attack is a key effectiveness and safety factor. The safe 
angles of attack and sideslip attainable in level and maneuvering flight 
often represent the important differences in both offensive/defensive 
capability and accident histories of otherwise equivalent aircraft. 1n 
the past generation of fighter aircraft, surprise loss of control depar¬ 
ture generally has been of concern as a key maneuver-limiting factor. 
The more recent generation has emphsized design for departure resistance 
but not necessarily prevention. Many operational pilots desire the 
capability to depart the aircraft on command as a last-ditch defensive 
maneuver, since departure is generally violent with unpredictable 
results from the attacker's viewpoint. The key element then revolves 
around departure resistance (or susceptibility), warning, severity, and 
recoverability. Accordingly, the issue of improved utility and safety 
at high angles of attack has been re-recognized as an area of fruitful 
research and development. 

The recent upsurge in interest has resulted in numerous studies 
devoted to identifying the safety-related maneuver-limiting phenomena 
inherent in: 

• Buffet 

• Pitchup or "dig-in" 

• Roll reversal 

• Wing rock 

• Nose slice 

• Rolling departure — snap roll 
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Some of these phenomena are perturbations about a steady flight 
condition (e.g., buffet, wing rock), which not only limit tracking 
capability (as shown in Reference 1) but also may serve as warning of 
approaching “departure" from controlled flight (e.g., nose slice, roll¬ 
ing departure). The more interesting and mysterious of these phenomena 
are the latter, and considerable effort has been expended in correlating 
the angles of attack at which such behavior occurs with various aero¬ 
dynamic stability and control parameters (e.g., References 2-4). Such 
open-loop "correlations," while useful, are not so satisfying or in¬ 
structive as the more positive identification of causal factors offered 
by the methods developed in Reference 5. 

The difficulty is in identifying causal relationships in the complex 
Interactive pilot/vehicle situation during uncoordinated flight, where 
aerodynamic cross-coupling can compound vehicle dynamic characteris¬ 
tics. The previous (Reference S) analysis and pilot simulation showed 
the static aerodynamic cross-coupling coefficients and N£ to contri¬ 
bute to closed-loop nose-slice departure susceptibility and severity in 
the A-7 aircraft. It represented an Initial successful and promising 
new attack on an old problem. But configuration differences among air¬ 
craft types, models and even loadings have been observed to produce 
gross differences in behavior through changed aerodynamics and inertial 
characteristics. Thus it is desirable that the same methods be applied 
to additional high-performance fighter aircraft having widely differing 
high AOA handling characteristics to see if further cause-effect rela¬ 
tionships can be identified. 

Accordingly, the stated or implied goals of this program are to: 

• Identify key design parameters that limit hlgh- 
angle-of-attack maneuverability for contemporary 
high-performance attack and fighter-type aircraft. 

• Postulate fundamental aerodynamic and control 
system design methodologies that will alleviate 
the limiting conditions. 

• Formulate handling qualities requirements for high 
AOA maneuvering flight to be incorporated in MIL- 
F-8785C, the military flying qualities specifica¬ 
tion for piloted airplanes. 
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B. TE CUBICAL APPROACH 


The research encompassed four major technical areas. The first was 
devoted to development and validation of aerodynamic models for two 
fighter aircraft having significantly different high AOA maneuver- 
limiting characteristics. The two aircraft selected, the F-4J and the 
F-14A, were found to represent almost opposite extremes in their range 
of departure susceptibility and severity. Attention was then turned to 
investigation of maneuver-limiting characteristics of each airframe and 
identification of causal relationships based on quasi-linear* analysis 
at symmetric and asymmetric ((3 i 4 0) flight conditions. The causal 
parameters were varied to alter the dynamic characteristics of each 
airframe at high AOA. Predicted characteristics were checked utilizing 
complete six-degree-of-freedom <6 DOF) models with nonlinear 
aerodynamics, and any differences between the quasi-linear frozen point 
analytic predictions and the nonlinear model results were resolved. 

The second task involved development of methods, criteria, and an 
associated pilot rating scale, for evaluation of handling qualities when 
approaching controllability limits. Previous analysis and simulation 
(Reference 5) has demonstrated that the Cooper-Harper handling quality 
rating scale is inapplicable for departure/recovery flight situations. 
Therefore, a new pilot rating scale is required specifically for high 
angles of attack. 

The third teak, pilot simulation, was performed to evaluate the 
influence of intentionally varied maneuver-limiting characteristics 
(wing rock, nose slice, etc.) under different normal flying situations 
such as training flights and air combat tracking tasks; to identify key 


A 

It is recognized that these linear analysis techniques are applied 
to phenomena which may be nonlinear in nature. However, the intent of 
the study is to investigate in an understandable manner the conditions 
which can lead to or precipitate rapid changes in motion, as opposed to 
analyzing the fully developed large angle phenomena. Therefore, piece- 
wise linearity is appropriate — but care should be taken in attempting 
to extrapolate the results. 
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flying quality parameters in terms of departure resistance, warning, 
severity, and recovery for possible Inclusion in the flying qualities 
specification; and to exercise and refine as necessary the new pilot 
rating scale. 

The fourth task involved assessment and correlation of results of 
the simulation with other applicable results and data; and formulation 
of high AOA maneuvering flight generalized design guides and flying 
qualities criteria for incorporation into MIL-F-8785C. 







SECTION II 


AIRCRAFT HIGH AOA DTMAKIC CHARACTERISTICS 

Thla section presents condensed results of an extensive analysis of 
the high AOA dynamic characteristics of the two example aircraft. The 
F-4J is treated first, followed by the F-14A. The high AOA character¬ 
istics of each aircraft are described as determined in flight test and 
operational deployment. The analytical models are briefly summarized 
and example comparisons shown to demonstrate the match obtained with the 
actual aircraft. 

Results of linear, frozen-point dynamic analysis at a series of 
symmetric and asymmetric (B 4 0) flight conditions are summarized. This 
analysis provides insight to potential causal factors behind undesirable 
characteristics at high AOA, showing examples of advantages and limita¬ 
tions of such analysis whan applied to situations in which aerodynamics 
are obviously quite nonlinear. In particular, this section Identified 
lateral-directional static aerodynamic cross-coupling as a key factor in 
determining high AOA stability and controllability characteristics. 

A. F-4J ANALYSIS 

The F-4J model of the F-4 family was selected because a previously 
validated aerodynamic model (Reference 6) was in hand. However, all 
versions of the aircraft with an unslatted leading edge wing are 
reported to exhibit quite similar high AOA flying characteristics: 
increasing AOA produces buffet of increasing severity, increasingly 
adverse yaw, wing rock, mild pitch-up, stall, and finally departure. A 
representative trimmed lift curve with onset of these handling phenomena 
is shown in Figure 1 from Reference 7. Above about 11 deg AOA, aileron 
adverse yaw becomes significant and rudder is used for roll control. 

Wing rock appears primarily as a rolling-sideslipping motion. It 
is often described by pilots as a divergent dutch roll which is easily 
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MODEL F-4 AIRCRAFT 
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Figure 1. Basic F-4 Trimmed Lift Curve (from Ref. 7) 

aggravated or caused by pilot lateral stick inputs. The angle of attack 
for onset and the severity actually depend upon several facto.. ,, includ¬ 
ing c.g. location, roll inertia or loading (clean vs. external stores), 
gear and flap settings, etc. It is generally more pronounced at light 
weight, with low roll inertia, and at forward c.g. The oscillations 
generally require 10 sec or more of sustained high AOA to develop. 

Amplitudes can be large enough (±30 deg $ and ±10 deg ?) to interfere 
with tracking performance, or to be highly uncomfortable when escorting 
considerably slower aircraft. 

The mild pitch-up is alleviated by a pitch rate damper but does 

result in some lightening of stick forces. It has been described by 
some pilots as a "dig-in" during slow turns. Possibly it is more of a 
stall warning cue than a maneuver-limiting problem. 

Departure is described variously as nose slice, roll departure, or 
both. Its onset also appears to depend upon c.g. location and aircraft 

inertia variations. Although shown in Figure 1 to occur at considerably 

higher AOA than wing rock, departure is not necessarily preceded by wing 
rock — especially if AOA is rapidly increased. 
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1. Analytic Model 


The F-4J aerodynamic model (Reference 6) was derived from three 
excellent but separate data sources (References 8-10). It was assembled 
to model as closely as possible the F-4J sub- and transonlcally for a 
moving-base simulation used to train U.S. Navy pilots In air combat 
maneuvering. It therefore had to be a very good representation of the 
complete aircraft in high AOA maneuvering up through stall. The instruc¬ 
tor pilots pronounced the simulation "credible and realistic of actual 
F-4J handling characteristics." 

The history of development of the complete aerodynamic model and the 
data package is presented in Part III, Appendix I. In brief, the aero¬ 
dynamic coefficients represent the following flight regime: 

0 < a < +110 deg 
-30 < 8 < +30 deg 

For this analysis the look-up table data were restricted to low Mach 
(<0.4) and a single altitude (h • 15,000). Weight and inertia charac¬ 
teristics representative of partially full internal fuel tanks, an empty 
600 gal centerline tank, and missile pylons on wing stations 2 and 8 
have been asaumed. 

Conventional 6 DOF equations plus kinematic terms are employed with 
the moment equations in body centerline axes, the force equations in 
wind axes, and aircraft orientation angles in standard Euler axes. The 
nonlinear equations of motion, aerodynamic force and moment equations, 
auxiliary equations, etc., are detailed in Appendix I (Tart III). 

A block diagram of the basic flight control and augmentation mech¬ 
anization is presented in Figure 2. Aircraft control is exerted through 
horizontal e'cabllator ( 5 8ta b)» rudder (6 r ), aileron (6 fl ), and spoiler 
(5. n ) surfaces. The aileron and spoiler systems are interconnected 
Buch that lateral deflection of ths control stick produces downward 






a) Longitudinal 




Figure 2* F-4 j Flight Control System 






deflection of one aileron and upward deflection of the opposite spoiler. 
Stability augmentation is provided in all three axes, but it Is common 
practice for the pilot to turn the roll damper off during high AOA 
maneuvering. The complex artificial feel system was not modeled in this 
analysis and simulation. 

2. Model Validation 

As a part of the checkout and acceptance test, the simulation based 
on the Reference 6 aerodynamic data was flown through various offensive 
and defensive air combat maneuvers, stalls, and departures by Navy in¬ 
structor pilots, who indicated that it adequately represented the F-4J 
handling and performance. This provided the first gross validation of 
the aerodynamic model. 

Actual flight test traces of F-4J high AOA dynamic responses are not 
available to compare against those of the mathematical models. The 
available high AOA flight traces are from an F-4E stall/post-stall 
flight test (Reference 11). However, Reference 12 indicates that all 
hard-wing models of the F-4 have approximately the same stall/departure 
characteristics. The F-4E has a nose section approximately 5 ft longer 
than the F-4J; this could cause some difference in side force due to 
asymmetric vortex shedding at high AOA (Reference 13). In addition the 
flight test vehicle was equipped with a spin chute, had a reinforced aft 
fuselage structure to handle the spin chute loads, and is presumed to 
have offsetting ballast in the nose. Thus, its pitch and yaw Inertia 
characteristics are quite different from the F-4J, as shown in Table 1. 
This difference is assumed to have minor Influence on the sequence of 
stall/departure characteristics with increasing AOA. The larger yaw 
inertia and longer nose (greater vortex shedding) of the F-4E could be 
offsetting departure susceptibility factors. But the inertia differ¬ 
ences should produce measurable differences in dynamic response parame¬ 
ters such as dutch roll frequency, response to control Inputs, and 
angular rates associated with departure. For example, scaling parame¬ 
ters by the ratio of inertias gives ratios of ti> ap and utj of 1.11. 
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TABLE 1. COMPARISON OF F-4E AND F-4J 
WEIGHTS AND INERTIAS 


F-4E (Reference II) F-4J 


<£. stores 

Empty tank 

Empty tank 

Wing stores 

Pylons 1, 2, 8, 9 

Pylons 2, I 

W (lb) 

40,000 

37,000 

c.g. (X c) 

28.1 

29.3 

I x (slug-ft 2 ) 

27,500 

23,850 

I y (8lUg-ft 2 ) 

157,000 

127,400 

I z (slug-ft 2 ) 

180,600 

146,000 

I X2 (slug-ft 2 ) 

5,500 

2,210 


Example comparisons between the F-4E test results and the nonlinear 
6 DOF F-4J math model are shown In Figures 3-5. Figure 3 shows the F-4E 
at a nearly steady AOA averaging about 23 deg. The rudder trace indi¬ 
cates that the yaw damper is on. The wing rock starts as a divergent 
dutch roll but then appears to limit at approximately 40-50 deg roll 
amplitude with a $/8 ratio between 2 and 3. Also the divergence grows 
linearly, instead of exponentially as it would at constant c • During 
the large, constant-amplitude 4> and 8 oscillation, the AOA trace shows a 
small oscillation at twice the frequency of the lateral oscillation. 

Figure 4 is the F-4J simulation trimmed to a Q ■ 23 deg and disturbed 
by a lateral stick pulse. Here the dutch roll also diverges almost 
linearly for about three cycles and then becomes limited in amplitude. 
The $/& ratio is about 3 during the limit cycle. An oscillation devel¬ 
ops in the AOA trace, again with a frequency double that of the lateral 
oscillation. The characteristics are remarkably similar except for the 
period of the oscillations. The higher frequency of the F-4E oscilla¬ 
tions is most likely due to a higher dynamic pressure. In Figure 3 the 
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Figure 4, F-4J Model Wing Rock at a 0 - 23 deg, 6 
Fitch and Yaw SAS On 









F-4E FLIGHT TEST 



Figure 5. Departure During Left Turn 


aircraft was pulled up to a * 23 deg within about 10 sec and then held 
at this AOA for about 13 sec. The initial Mach number was not indi¬ 
cated. In Figure 4 the model was trimmed at an initial speed corre¬ 
sponding to a Q - 23 deg. If all other factors were exactly the same, 
the difference in frequency shown here would be accounted for with only 
a 30 percent difference in speed. 

The computer-generated images of Figure 5 show a comparison of the 

F-AE aircraft and the F-4J model departures from a wind-up turn to the 

left. The departure onset is quite similar. Both depart with a nose 
slice away from the turn accompanied by a roll oscillation. The F-4E 
ends up with more of a rolling departure, while the F-4J ends with more 
of a yaw departure (nose slice). This difference in final motion is 
probably influenced greatly by the differences in inertias. The F-4E 
has 24 percent higher yaw inertia and 7.3 percent higher ratio of yaw to 

roll inertia. Thus it could be expected to be more resistant to yaw 

departure and prone to roll departure. 
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Since our main interest is in modeling and establishing causal rela¬ 
tionships leading to high AOA maneuver-limiting factors rather than 
modeling post-stall gyration and spin, the foregoing match of the non¬ 
linear wing rock behavior and nose slice onset is considered to be one 
validation of the aerodynamic model. Appendix I of Part III and por¬ 
tions of the analysis and simulation to be presented later offer addi¬ 
tional validation. 

One key aspect is that the aerodynamic model demonstrates both limit 
cycle wing rock and nose slice departure without introducing artificial 
hysteresis in aerodynamic coefficients (Reference 14), drastic changes 
in roll damping (Reference 10), or contrived elevator control asymmetry 
(Reference 15). 

3. Dynamic Analysis 

In high AOA maneuvering flight it is not uncommon for the aircraft 
to be in asymmetric (6 j 0) flight — either intentionally or uninten¬ 
tionally. Previous studies (e.g., References 5, 16-18) have shown that 
static aerodynamic cross-coupling due to sideslip can have significant 
influence on aircraft stability and dynamic characteristics and possibly 
the departure characteristics. For example, Figures 6-8 reflect the 
pitch, roll, and yaw aerodynamic moments, respectively, for our F-4J 
model as a function of a and 8. For simplicity the aerodynamic moments 
were assumed symmetric with sideslip. Over the region 15 < a < 25 deg 
pitching moment is significantly Influenced by sideslip, and rolling and 
yawing moments due to sideslip vary greatly with AOA. This is the same 
AOA region in which the various maneuver-limiting characteristics of the 
F-4 are exhibited. Thus, in the analysis to follow, attention will be 
devoted to examining the possible influence of static aerodynamic cross- 
coupling on high AOA maneuver-limiting factors. Aircraft open-loop 
stability characteristics will be analyzed first, and then cloBed-loop 
piloted control aspects will be determined. 
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Figure 8. F-4J Yew Coefficient Variation 
with a, 8 
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a. Open-Loop Dynamic Characteristics 

Figure 9 preaenta a locus of longitudinal short-period an 1 lateral 
dutch roll root migrations with AOA in wings-level, symmetric, 1-g 
flight at an altitude of 15,000 ft. These dynamic parameters represent 
symmetric, linearised, fixed-operating-point conditions. With 

increasing AOA the dutch roll mode is seen to progress from a lightly 
damped stable mode to oscillatory divergence and, finally, to a pair of 
aperiodic divergences while the longitudinal short period remains 
relatively unchanged. 

Also shown in the figure is the AOA at which the dynamic stability 
parameter Cng^yn becomes zero. At lower AOA the parameter is positive, 
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- 1.0 -.5 0 ~ ,5 ' 1.0 

<7 

Figure 9. F-4J Lateral-Longitudinal Short-Period Root 
Migration with a; 0 O “ 0 

Indicating directional dynamic stability. At high AOA the parameter is 
negative, indicating directional dynamic instability. The AOA regions 
below and above Cnpjy n - 0 are identified as wing-rock and nose-slice 
regions, which correlate quite well with the AOA rangeB Identified with 
these characteristics in the F-4E flight test, Figure 10 (from Refer¬ 
ence 19). 

The influence of non-zero sideslip at trim on the dutch roll and 
short-period frequencies and damping 1 b reflected in Figure 11. The 
6 DOF aircraft was trimmed to l.S and 5.5 deg sideslip and linearized 
transfer function parameters obtained. These indicate that the dutch 
roll is destabilized and the short-period stabilized by steady sideslip. 
Thus there is an apparent Interchange of damping between the two moden 
while the overall system damping remains essentially constant. 

Unfortunately, comparison of dynamic characteristics predicted by 
linear frozen-point analysis in Figures 9 and 11 with the actual flight 
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Aim Norma! 


Symbol 

Gross Weight 
(lb) 

Load Factor 

(g) 

O 

40,500 

1,0 

□ 

40,000 

1.0 

A 

39,400 

1.0 

T 

39,200 

1.0 

♦ 

39,000 

1.5 

1 

39,200 

2.0 

tv 

37,100 

1.0 

A 

36,500 

1.0 

k 

36,200 

1.0 


35,700 

1.0 


Aim Normal 


Symbol 

Gross Weight 
(lb) 

Load Factor 

( g ) 

★ 

35,200 

2.0 

%r 

34,800 

3.0 

□r 

40,000 

1,0 

A 

39,700 

2.0 

V ' 

38,700 

2.0 


43,500 

1.0 

A 

43,200 

1.0 


42,700 

2.0 


42,200 

3.0 


41,700 

3.0 

0 ' 

41,300 

1.0 


a(deg) 



26 - 

a at Onset 
of Nose Slice 

21 - 


16 - 


26 - 

a at Heavy 
Wing Rock 

21 - 

(±20deg roll) 

16 - 



0 

9 
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a at Onset 
of Wing Rock 
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22 
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A 
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Figure 10. F-4E Stall Approach Characteristics 
(From Ref. 19) 
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Figure 11. F-4J 6 DOF Linearized Equation; Lateral- 
Longitudinal Root Migration with a and 6 


traces of Figure 3 or nonlinear 6 DOF model traces of Figure 4 shows 
several significant differences. The wing rock linear analysis predicts 
a dutch roll divergence of increasing amplitude as a, 3, or both are 
increased, and a longitudinal short-period with a frequency quite close 
to the dutch roll frequency, but heavily damped. The traceB of 
Figures 3 and 4 show the dutch roll divergence bounded in amplitude and 
a low-damped longitudinal oscillation at about twice the dutch roll 
frequency. 

Figure 12 shows a similar set of nonlinear time response character¬ 
istics for the open-loop 6 DOF nonlinear model trimmed to a Q - 21 deg 
and 3 Q ■ 0 deg and excited with an aileron pulse. From Figure 12 it is 
apparent that the longitudinal oscillation results from rectifying the 
lateral oscillation and the sharp peaks of pitching acceleration 
coincide with sideslip passing through zero. Thus, the nonlinear-type 
behavior appears to be caused by lateral-longitudinal coupling due to 
sideslip. 

A nine-by-nine matrix of aerodynamic and kinematic terms for 
coupled, non-symmetric flight obtained from the partial derivative 
expansion of the nine equations of motion is presented in Figure 13. 
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The major derivatives and kinematic terms are Identified in literal form 
and, for comparison with Figure 12, have been evaluated at a Q « 21 deg 
and 3 0 ■ 1.5 deg. The principal cross-coupling at this small sideslip 
angle Is seen to be due to Mg, «C £ and Nq. The Influences of these 

derivatives on actual vehicle dynamic response can best be seen through 

the use of eigenvectors and vector polygons. 

Figure 14 contains the eigenvectors for the dutch roll and short- 
period modes. These show both modes to be dominated by roll, side¬ 
slip, and pitch, with the difference between the modes being primarily 

motion phasing and damping. For example, in the dutch roll, q and 8 are 
nearly 180 deg out of phase, while in the short-period they are nearly 
in phase. However, the short-period is so well damped (C ap * 0.6) that 
It would Influence only the first second or two of the time responses 
shown in Figure 12. Therefore, the observed time responses must be due 
entirely to the dutch roll. 

Vector polygons for the four principal motion equations are pre¬ 
sented in Figure 15. Each polygon represents a force or moment equa¬ 
tion. The polygon's sides show the contributions of each motion varia¬ 

ble to the total acceleration. Vector lengths are a product of the 
appropriate stability derivatives and eigenvectors evaluated at the 

dutch roll frequency. The predicted motion is mostly rolling (p), with 
• « 

sideslip (@), and pitching (q) about the same magnitude, and with com¬ 
paratively little yawing (r). The polygons for a and u are so small 

that they can be neglected. The vector polygons provide the first 

insight to the caused for the change in modal response characteris¬ 
tics. The p vector polygon shows phasing of the ta , 0 vector opposite 
to that for the <£pp vector. ThuB the aerodynamic cross-coupling 
coefficient, in effect opposes the normal aerodynamic damping 

term, tp, and therefore increases the rolling tendency of the aircraft. 
Since the aerodynamic cross-coupling coefficients Increase in magnitude 
as the trim sideslip increases, one can readily visualize that at higher 
sideslip angles, to. can actually become greater than £p. Then one 
could expect the model to exhibit a roll divergence for large a and 3 
perturbations. 
















p Equation 


p = X.0 jfi + t'a a + <£rf + £pP 



Figure 15. F-4J Force Vector Folygons; 

«o “ 21 d «8. 6 0 “ 1,5 deg 

The 3 vector polygon Indicates that this vehicle motion derives pri¬ 
marily from the kinematic coupling of roll rate about the body center- 
line axle at the elevated AOA. 

In the r vector polygon the p vector la missing because the aero¬ 
dynamic coefficient Np is negligible at this flight condition. The 
vector N£a has a relatively large component in phase with N£r so that 
the aerodynamic cross-coupling from longitudinal into lateral augments 
the natural yaw damping and thus contributes to reducing body-axis 
yawing motion. All of this contributes to the dominance of roll and 
sideslip in the Figure 12 wing rock. 

Finally the q vector polygon shows that the MgB vector is of oppo¬ 
site phase and is larger In magnitude than the damping sector derived 
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from the aerodynamic derivatives Mq + M<J. Thus the sideslip-induced 
pitching moment overpowers the normal aerodynamic damping and forces the 
double frequency oscillation observed in the time traces. Mg changes 
sign with 0 (Figure 6). Since pitch inertia is considerably larger than 
roll inertia for this aircraft, a transfer of energy from the unstable 
lateral mode into the stable longitudinal mode may be a contributing 
factor in bounding the lateral-directional mode divergence. The influ¬ 
ence of reduci. 1 ; is shown in Figure 16. In this case, is small 
but not zero. Comparing these traces with those of Figure 12 demon¬ 
strates the strong influence Mg has in bounding the lateral divergence. 

It is concluded for our model that the roll/sideslip oscillation 
starts as a divergent dutch roll with small sideslip excursions. As 
sideslip ‘increases, coupling causes a nonlinear behavior which bounds 
the divergence. Thus F-4J wing rock, per se, 1 b a nonlinear phenomenon 
dominated by static aerodynamic lateral-longitudinal cross-coupling. 
Further, the vector polygons show that except for the damping 
derivatives do not have a significant influence on high AOA dynamic 
characteristics. 

This analysis has also shown that one must exercise caution in the 
use or Interpretation of open-loop transfer function parameters even 
when obtained from 6 DOF fixed operating point conditions. Such param¬ 
eters are limited to relatively small sideslip conditions such as onset 
to wing rock or departure and hence better reflect susceptibility to 
these maneuver-limiting factors rather than describing the resulting 
vehicle motion. 

b. Closed-Loop Dynamic Characteristics 

The cloBed-loop dynamic characteristics of any system are strongly 
Influenced by transfer function numerator roots (zeros) in or near the 
range of the desired bandwidth of control. When controlling any motion 
variable with any controller, the numerator zeros of that transfer func¬ 
tion attract denominator (open-loop) poles. Hence these zeros alter 
closed-loop stability and modal response characteristics of the vehicle 
in proportion to the tightness (gain) of the loop closure. Thus certain 
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zeros can become high AOA maneuver-limiting factors in that they cause 
or influence closed-loop roll reversal, wing rock (lateral PIO), and 
departure (nose slice or roll). This observation applies to loops 
closed both by the pilot and by automatic feedback systems. The key 
numerators for piloted control are N t r » an< * ^5 s tab* T ^ e 

augmentation systems are of little concern due to limited SAS authori¬ 
ties and very low loop gains at high AOA and low dynamic pressure. 


Lateral stick control of roll . Under normal low AOA flight 
situations, the 3 DOF equations produce a transfer function 


$(s) 



N ^stk " + + u $] 


A - (s + 1/T g )(s + 1/T r )(b2 + 2c d o) d s + to d ) 


with values of and 5^ in the vicinity of those for the dutch roll 
mode, and C d ^ d . With this loop closure the zeros attract (modify) 
the dutch roll, causing the closed-loop root locus (with increasing 
feedback gain K) to take on a form typified by the sketch below. 
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2 

As AOA Is increased, generally decreases and can become negative 
in value. When this occurs, form 

“U ■ A +(* + ^r)(* + t^) 

where one root is positive and one le negative. With the high AOA shift 
in open-loop denominator roots, the root locus often takes on the form 
shown below. 



In this case the spiral root, -1/T g , is driven toward the right 
half-plane (RHP) zero and a first-order instability (divergence) re¬ 
sults. The rate of divergence depends upon how far the zero lies in the 
RHP and how tightly the loop is closed. 

Migration of the F-4J open-loop zeros with AOA at zero sideslip 

is reflected by the solid lines in Figure 17, which indicate an almost 
2 

linear decrease in “a (at constant damping) with increase in AOA up to 
• 2 

a - 18 deg. Above this angle, is negative and the RHP root location 
is quite AOA-sensitive. Also noted on the figure are the regions asso¬ 
ciated with wing rock and nose slice from Figure 10. Wing rock occurs 
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2 

in the region where is small; nose slice occurs in the region where 
1 /T^ is large and negative. 

2 

Actually, negative ^ reflects steady-state roll opposite to that 
commanded by stick deflection. This is generally due to large adverse 
aileron yaw, since (after Reference 20) 



where + 11.66 ep . At high AOA, Ng is small, <£3 is large 

negative, and N$ ^ determines the sign of Adverse (negative) Ng stk 
then tends to produce negative and the RHP zero. In physical terms, 
adverse yaw produces sideslip and, in turn, roll reversal due to 
It thus becomes apparent that the nose slice region identified in Fig¬ 
ure 17 may be related to strong adverse aileron yaw and the magnitude of 
the RHP zero. 

2 

The dashed lines of Figure 17 show the values of ^ and 2C^ta^ ( 0 r 
1/T^, l/T,^) for e o - 5.5 deg. Whereas increasing AOA at 8 • 0 causes 
u>| to decrease steadily and finally become negative, when 8^0 the 
aerodynamic cross-coupling causes to stay positive nnd close to the 
dutch roll mode but the damping ratio, becomes negative (see sketch 
below). 


Lateral - Directional 
Root Locations for: 
a 0 = 21° 

/3 0 = 5.5° 


W SH 




u.'<j 
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The promixity of and means low residue for the dutch roll mode 
In roll control with stick when 0^0. Thus sideslip should reduce the 
roll content of the unstable dutch roll mode; this Is consistent with 
the observed bounding of the roll oscillation divergence in Figures 4 
and 16. 

Rudder control of roll . When the region of adverse aileron yaw Is 
reached, roll control is maintained via rudder. This produces yaw and 
sideslip-induced roll in the same direction and thus much more rapid 
turning than just yawing moment and side force would produce. The 
rudder-to-roll numerator, , is of the form 



One root Is always positive (RHP) as shown in Figure 18. The RHP zero 
results from adverse roll due to It is apparent from Figure 18 
that the zeros of are not greatly influenced by either AOA or side¬ 
slip. 

Again this control structure portends closed-loop instability since 
some pole must be driven toward the RHP zero. However, these zeros lie 
at high frequencies relative to the low-frequency bandwidth normally 
achieved with rudder control via the pilot's legs. Rudder-to-roll con¬ 
trol handling qualities are normally dominated by the low-frequency 
spiral and roll subsidence modes, and at high AOA these remain stable 
and well behaved. While this type of control is "unnatural," it does 
not result in roll stability problems at low gain closures. 









Figure 18. Root Migration with AOA{ F-4J 


Stabllator control of pitch . Closed-loop control of pitch attitude 

becomes pertinent to maneuver-limiting If and when the numerator for 

0 

stabllator control of pitch, N S stat) » h*s aeros In the RHP. At aero 
sideslip the transfer function is of the form 


0(s) 


stab 

Z 




A ■ [s^ + 2CpWpS + <*) p ] [ s + 2i; B p(ij B p8 + t*>§p] 
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On the F-4J, I/Tq^ decreases with increasing AOA and finally couples 
with 1/Tg^ to form a complex pair which almost cancels the longitudinal 
phugoid mode. Thus these zeros remain well behaved. However, when 
there is sideslip so that static aerodynamic cross-coupling becomes 
significant, the numerator contains additional zeros related to lateral- 
directional modes. 

It was discovered (Reference 5) that such coupling in the A-7 air¬ 
craft produced an RHP zero, I/Tq^, and subsequent pilot control of pitch 
attitude via the stabilator caused a closed-loop directional divergence 
which has the appearance of a nose slice. A similar situation exists 
in the F-4J aircraft. Figure 19 presents pole and zero locations for 
a 0 - 23 deg and $ 0 » 5.S deg. The coupled 6 DOF equations (Figure 19a) 
produce two real zeros (one RHP, one LHP), one complex pair which nearly 
cancels the phugoid, and one higher-frequency complex pair. 

The source of the low-frequency zeros can be determined by elimin¬ 
ating the velocity equation, which should reduce the phugoid mode to 
a first-order pole at the origin and eliminate the speed numerator zero, 
l/Te^. It may be observed in Figure 19b that the only change is in the 
phugoid dipole pair which became first-order as expected. 

Next, as in Reference S, the cross-coupling terms and N£ are set 
to zero in the 6 DOF equations to identify the lateral-directional 
modes. This is shown in Figure 19c where the pole-zero cancellations 
show the two real zeros to be related to the roll subsidence, 1/T r , and 
spiral, l/T # , modes and the higher-frequeucy complex zeros to be related 
to the dutch roll. 

Figure 20 is a system survey for closure of the pitch attitude loop 
for the augmented (SAS on) airframe at this same trim condition. The 
lower figure is a root Iocub and the upper figure is a Bode-siggle 
(Reference 20). The root locus shows the migration of roots from the 

'''The slight difference In airframe dynamic parameters between 
Figures 20 and 11 indicates the very minor Influence of the SASs at high 
AOA and low dynamic pressure conditions. 
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open-loop poles toward the open-loop zeros as gain Is increased. It 
provides information concerning closed-loop frequency and damping as a 
function of gain but does not provide Information concerning sensitivity 
to gain change. The Bode-siggie presents the open-loop transfer func¬ 
tion amplitude and phase as functions of frequency where the amplitude 
is evaluated for s * jw (Bode) and s ■ o (siggle). The symbols iden¬ 
tify migration of complex LHP closed-loop pole frequency with increasing 
gain (Bode-root locus, Reference 20) and the 0 symbols similar migra¬ 
tion of closed-loop pules in the RHP. Closed-loop first-order roots 
lie at the intersection of the gain closure lines and the o amplitude 
curves. Thus these survey plots allow identification of closed-loop 
roots at a specific gain ilue, and sensitivity of root value to changes 
in gain. 

In order to provide any pitch attitude control whatsoever, the loop 
gain must he sufficiently high that the closure line crosses the Bode 
plot below the -20 dB low-frequency asymptote, i.e., with a gain 20 dB 
or greater. Generally the crossover will be at or above the short- 

period frequency — again in the vicinity of -30 to -40 dB. The root 
locus plot has tick marks showing closed-loop root locations for gains 
of 20 and 32 dB. This shove that the root moving toward l/Tg^ has 
already reached the RHP, and divergence, with the gain at 20 dB. By 
32 dB the first-order divergence has a time constant of 4 sec. Thus any 
attempted control of pitch attitude results in a divergence. This 
divergence la due to lateral-longitudinal coupling and is evidenced in a 
lateral motion. 

0 

Figure 21 shows values for the real zeros 1/Tg.j and 1/T R at @ 0 - 

5.5 deg and a Q from 21 to 25 deg. These Indicate that at non-zero 

sideslip and A0A Increasing above 20 deg 1/Tg^ rapidly moves to negative 

values which portend closed-loop first-order divergence time constants 
as low as 1 sec. Again this coincides with the nose slice region 
observed in flight test. 
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Figure 21. N$ First-Order Root Migration* 
with AOA; F-4J 


Coupling numerator* . when the aircraft ia in asymmetric, 
sideslipping flight, any control action in the longitudinal axis 
influences the lateral-directional numerator and vice versa. This can 
be illustrated as follows: 


Alrfrome 








The effective transfer function for each control loop becomes 


"Lab . 8 *.t.b + V ^.t.blt1.' 

A- -£- 

. *t..k + 

- -1- 

where N °stab*.tk is the numerator term due to cross-coupling effects.* 
Both the 0 and $ numerators are modified by the same coupling numerator; 
the extent of the influence depends upon how tightly the other loop is 
being closed. This Influence can be determined by treating the 
effective numerator as a closed-loop system and identifying the 
migration of numerator roots as the loop gain (Y^ or Y e ) is varied, 

e.g., the roots of the 6 numerator ars obtained from 

K l.Ub “ "tab + " ° 

1 + Y *%talXtk . 0 

N «stab 

The uncoupled numerator roots become the poles of this system and 
the coupling numerator roots are the zeros. The root migration is 

controlled by the pilot gain, Y^, For the Figure 20 example case of 
a Q m 23 deg, 0 O ■ 5.5 deg, the transfer function is 



t| 

The reader is referred to page 163 et seq. of Reference 20 for a 
thorough development and treatiae of multiloop control systems analysis. 
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where first-order factors (s + a) are represented by (a) and second- 
order factors [s + 2cto n s + u£] by [c; u n ]. Assuming Y^ to be a simple 
positive gain, this produces a root locus of the form shown in the 
sketch below. Thus pilot control of the roll attitude loop via aileron 
and spoiler causes the RHP zero of the pitch numerator to move further 
into the RHP because the root of the coupling numerator lies beyond 
1/T 03 . 



Similarly, pilot control of pitch attitude via Btabilator alters the 
roll numerator as 



stk 


n; 


stk 


-.0029(-1.35)(1.69)[.993; .102]Y e 

7M3TI: re7y • iwr t tgi t 


which, for Yg assumed a simple negative gain, produces a root locus of 
the form sketched below. This results in the RHP complex pair being 
driven to the real axis with one root eventually moving further into the 
RHP while the other moves into the LHP. 







Thus in each instance the influence of the coupling numerator is 
definitely not beneficial; the coupling can further degrade closed-loop 
control in both axes. 

Since F-4 pilots are normally instructed to use rudder to roll at 
high AO A, the Influence of this control technique may be observed from 
the transfer functions (again at a Q ■ 23 deg, 0 ■ 5.5 deg): 



Use of rudder to roll modifies the pitch numerator in much the same 
manner as doeB use of aileron and spoiler. This can be aeen by 
comparing values of the coupling numerator terms. On the other hand, 
pitch control has little Influence on the rudder roll numerator because 
nearly all numerator and denominator roots cancel. Thus it appears that 
this control strategy produces the most consistent closed-loop lateral- 
directional dynamic properties at high AOA and non-zero sideslip, which 
also might be a reason pilots do not mind using this technique. 

Summary . This analysis has shown several areas in which control- 
numerator dynamic parameters can become maneuver-limiting factors for 
the F-4J. These limits are generally associated with zeros migrating to 
the RHP at high AOA in asymmetric (0 4 0) flight. 

The RHP zero of i* due to adverse "aileron" yaw nega¬ 
tive) and results in roll reversal. The attendant sideslip causes the 
4 zeros to shift toward the dutch roll poles. This proximity reduces 
the roll contribution in the dutch roll mode and produces a nonlinear 
bounding of dutch roll (or wing rock) oscillation amplitude. 

One RHP zero exists in N$ whether the aircraft is in asymmetric 
°r 

flight or not. However, that root is at a frequency sufficiently above 
the region in which rudder control is generally exerted that it does not 
pose a closed-loop stability problem. 
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An RHP zero can occur in N 6 sta ^ due to aerodynamic cross-coupling 
with non-zero sideslip* This zero results in a first-order divergence 
mode if stabilator control of pitch attitude le attempted. 

Asymmetric flight also produces lateral-longitudinal control cou¬ 
pling that further aggravates these undesirable RHP zero locations of 
the key € and $ numerators. 

Thus the aerodynamic cross-coupling associated with sideslip tends 
to dominate the F-4J high AOA flying qualities. Values of the cross- 
coupling derivatives evaluated at B 0 ■ 5.5 deg and a range of AOA are 
shown in Figure 22. Comparing these with the root migrations in 
Figures 9 and 11 demonstrates their significance. The regions where 
large shifts in dynamic parametars occur coincide with large magnitude 
of cross-coupling derivatives. Also, it should be noted that the 
magnitude of the cross-coupling derivatives is directly proportional to 
the magnitude of sideslip. 

The foregoing analytic results must be viewed with some caution 
since it was shown that linearized, fixed-operating-point dynamic 
parameters do not necessarily reflect the actual open-loop dynamics of 
the nonlinear airframe. On the other hand, the apparent correlation 
between bounaing of dutch roll divergence and movement of zeros with 



Figure 22. F-4J Static Aerodynamic Cross-Coupling 
Derivatives; 0 Q « 5.5 deg 
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sideslip, and Che AOA region for nose slice observed In actual flight 
with the region where the RHP zeros of 8 cab an< * N ^stk are *- ar 8® 18 
quite encouraging. 


B. F-14A ANALYSIS 

The F-14A configuration does not have wing leading-edge slats. This 
aircraft has variable wing sweep; at M < 0.5 the wings are in the full 
forward (A ■ 22 deg) position. The unaugmented airframe Is classified 
as departure resistant, but it can be departed (References 21 and 22). 
It remains laterally stable and longitudinally controllable at full aft 
stick (AOA approximately 50 deg with wings forward). It has a mild 
dutch roll Instability starting at about 15 deg AOA. The oscillations 
are aggravated if the pilot attempts to oppose the motion with either 
lateral stick or rudder. The dutch roll mode becomes stable again at 
about 23-25 deg AOA. Roll reversal due to adverse yaw from the aileron 
and differential horizontal stabilizer also starts at about 18 deg 
AOA. Prolonged deflection of the lateral stick at a > 2u deg will 
result in a departure, characterized by a rapid increase in "adverse" 
yaw rate, which may develop into a series of uncommanded rolls with the 
appearance of snap rolls. At low speed, departure can occur at AOA as 
low as 21 deg. The departure characteristics are a function of wing 
sweep, c.g. location, and Mach number. The airplane will spin and has a 
non-recoverable flat spin mode. 

1. Analytic Model 

The F-14A aerodynamic model was derived from several sources (Refer¬ 
ences 23-27). The data bane was taken from the NASA/Langley Differen¬ 
tial Maneuvering Simulator (DMS) report of Reference 23. This fixed- 
base piloted simulation was considered by Grumman test pilots to provide 
a valid representation of high AOA dynamic characteristics observed in 
early flight tests. However, time traces obtained from ther<n same aero¬ 
dynamic data did not match available flight test records. Therefore 
the date were modified on the basis of information obtfined from other 


4 ?. 








referencesi discussions with F-14 aerodynamleiats at Grumman, and 
finally, as necessary to match available flight test data. 

The resulting aerodynamic model, together with supportive documents'-' 
tlon for each aerodynamic coefficient, is detailed in Part III, Appendix 
II. The data caver the range: 

0 < a < 55 deg; -20 < 6 < +20 deg 


at low Mach (0.2). At this low speed the variable-sweep wing is 
nominally in its full forward (A ■ 22 deg) position, so wing sweep 
Influence is not Included. 

Based upon the F-4J analysis, the key high AOA aerodynamic moments 
are Cy,(a,g), and C m (a,@). The roll moment (Figure 2.3) is seen 

to remain essentially constant for a > 15 deg and 3 < 10 deg, as opposed 

to the F-4J'n significant decrease in Gfc(a,B) fora > 15 deg (Figure 7). 

The yaw moment (Figure 24) ia similar to the F-4J's (Figure 8) for 
a < 17 deg, but for a > 17 deg both the slopes with respect to a and the 
maximum values with respect to 3 are approximately half those of the 
F--4J. The combination of less negative C n ^ and larger negative Cig pro¬ 
duces positive values of the dynamic stability parameter 0ngdy n through¬ 
out the 55 deg AOA range. 

The pitch moment with sideslip (Figure 25) is also similar to that 

of the I-4J over the range 10 < a < 25 deg, where it is negative. How¬ 

ever, above a - 35 deg, becomes positive. 

Weight, c.g., and moments and product of inertia for the clean (no 
external stores) aircraft were provided by Grumman. These are shown in 
Table 2 with F-4J values for comparison. All inertias are In the 
fuselage centerline reference axis system. Aircraft dimensions are also 
presented for comparison. 




















TABLE 2 


COMPARISON OF F-I4A AND F-4J WEIGHTS, 
INERTIAS, AND DIMENSIONS 


Parameter 

F-14A 

F-4J 

W (lb) 

46,950 

37,000 

c.g. (X c) 

16.0 

29.3 

I x (slug-ft 2 ) 

58,950 

23,850 

I y (alug-ft 2 ) 

225,600 

127,400 

I z (slug-ft 2 ) 

285,000 

146,000 

I xz (slug-ft 2 ) 

-3030 

2210 

A (deg) 

22 

45 

CM 

4-» 

*W 

to 

565 

530 

b (ft) 

64.1 

38.67 

C (ft) 

9.8 

16.04 


The Inertia differences between the vehicles are large. Although 
the F-14A la some 27 percent heavier than the F-4J with comparable 
fuel loading, the inertias range from 77 to 147 percent higher and the 
product of inertia is of opposite sign. This has considerable influence 
on the ratios of aerodynamic to kinematic terms in the model and greatly 
influences coupling effects, as will be shown later. 

Aircraft control is exerted through a horizontal stabilizer (6 Bta b), 
twin rudders (6 r ), spoilers and differential horizontal (6 D ) sur¬ 
face deflections. Simplified block diagrams of the flight control and 
augmentation systems are shown in Figures 26-28. The longitudinal sys¬ 
tem contains a feel system (spring cartridge, variable gearing ratio, 
bobweights), a limited-authority pitch rate SAS, and a surface actuator. 
However, the analysis and simulation are based on a simple fixed stlck- 
to-Burface gain and pitch rate SAS. 
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Figure 28. F-14A Directional Control System 


The lateral system (Figure 27) has a roll rate command augmentation 
system (CAS) in parallel with a mechanical stick-surface link. At high 
AOA the spoilers are ineffective and all roll control is obtained 
through differential horizontal, <5 q. The roll rate CAS operates only 
through the differential horizontal, and the electrical path gains are 
varied with AOA such that the CAS command and feedback fades to zero 
when 20 < a < 31 deg and, simultaneously, the electric stick-surface 
path fades in to cancel all but ±2 deg of mechanical path surface 
command. Thua, for a > 31 deg there is no roll rate damping augmenta¬ 
tion and the pilot can only command ±2 deg of differential horizontal 
surface. 

The rudder system (Figure 28) contains a fixed gain mechanical link 
from pedals to surface actuators, a conventional yaw damper with yaw 
rate and lateral acceleration feedback, and a lateral stlck-to-rudder 
croasfeed (SRI). At a > 10 deg the yaw SAS is cut out and the SRI faded 
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in. The SRI signal, passed through a first-order lag filter (break 
frequency 8 rad/sec), reaches full gain at a > 20 deg. Above a * 20 deg 
full lateral stick deflection produces 19 deg of rudder deflection. 

Thus for the AOA region of Interest, a > 20 deg, the airframe is 
essentially unaugmented in the lateral-directional axes and the control 
stick effectively moves the stabilator for pitching and twin rudders for 
rolling. 

2. Model Validation 

F-14A flight test data covering the appropriate AOA region were 
provided by the aircraft manufacturer (Reference 28). Unfortunately, 
the data were extracted from tests of three aircraft having slightly 
different weights, inertias, flight control systems, and flight test 
instrumentation. Since the flight test goals involved flight control 
and leading edge flap system development rather than Identification of 
basic airframe dynamics, most runs were contaminated by SAS effects, 
pilot control inputs, or configuration Influences. The model validation 
therefore was based on general response wefcthirvg of a few short time 
"windows" of bare-airframa response or rel«'r' -aly uncontaminated traces, 
plus observation that the sequence of dynamic characteristic changes 
with increasing AOA described previously was indeed achieved. 

Figure 29 presents results of a static longitudinal validation check 
based upon a 1-g stall in which AOA is increased stesdlly from 10 deg to 
almost 40 deg with a ramp horizontal stabilizer input. The circled 
points reflect the excellent agreement for horizontal stabilizer inputs 
required to achieve the noted AOA trim condition at 35,000 ft altitude 
with the aerodynamic model. 

Figure 30 is a short time-slice of the divergent dutch rolj. at 
approximately 15 deg AOA. Motion is excited by a lateral stick doublet 
which, through the SRI, produces both differential horizontal and rudder 
doublets. Both the yaw and roll rate SASs are then turned off, and the 
subsequent oscillation represents the "bare" airframe in the lateral- 
directional sense. As with the F-4J, the envelope of the oscillation 
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Figure 29. F-14A Longitudinal Static Trim Comparison 


































seems to grow linearly rather than exponentially, which Indicates 
nonlinear damping. AIbo, the AOA trace shows an oscillation at twice 
the dutch roll frequency. The motion on the trace indicates the 

pitch SAS was on and responding to the longitudinal oscillation. 

Figure 31 presents time traces for the F-14 aerodynamic model 
trimmed at 15 deg AOA and excited by aileron and rudder doublets similar 
to Figure 30. This also shows an almost linear, rather than 
exponential, growth in dutch roll amplitude. The frequency and cycles 
to double amplitude are the same as in Figure 30. The small-amplitude 
AOA oscillation also is precisely twice the frequency of the dutch roll. 

Figure 32 presents another short time-slice in which AOA is 

maintained at 18 ± 2 deg. In this angle of attack region the SAS gain 

schedule shows that the yaw SAS should be off and the roll SAS on. 

However, comparison of the yaw rate and rudder position traces indicates 
that the rudder is correlated with yaw rate at approximately the 
magnitude and phase angle that would result from the yaw damper. Since 
the rudder pedal and lateral stick are both constant and very small, the 
aircraft response is essentially that of the controls-fixed augmented 
airframe. This dutch roll oscillation is essentially the same, 
indicating that the yaw SAS is quite ineffective at this AOA and 

airspeed. 

Other comparisons which can be made among Figures 30, 31, and 32 are 


summarised below: 

Parameter 

F-14A 

Figure 32 Figure 30 

Model 

Figure 31 

Period of oscillation (Bee) 

4.5 

4.5 

4.5 

Cycles to double 

2 

2 

2 

$/S (deg/deg) 

2 

* 

2 

p/r 

25 

** 

22.5 

a (deg) 

18 ± 2 

15 

15 


A A A 

3 trace not available. r trace contaminated. 
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Taken together, the foregoing is considered an acceptable match between 
actual and simulated airframe dynamics in the region of 15 deg AOA where 
the dutch roll divergence and wing rock are most pronounced. 

Figure 33 shows a pullup initiated from wings-level flight in which 
AOA is steadily Increased from 13 to 38 deg at about 0.75 deg/sec. A 
lateral stick oscillation is introduced at about 15 deg AOA. The stick 
is then centered and the SASs turned off at about 18 deg AOA. A small 
rudder ramp is then introduced, possibly to counter the Blow roll-off in 
bank angle. Simulation traces for a matching pullup with similar stick 
and rudder inputs are shown in Figure 34. Both Figures 33 and 34 show 
a divergent dutch roll. However, the amplitude appears to increase 
linearly rather than exponentially, again reflecting nonlinear damping. 
Both sets of traces show the period of oscillation to increase from 
about 4.5 sec to over 5.5 sec as AOA approaches 30 deg. Again there is 
good agreement between the actual flight and simulation responses. 

3. Dynamic Analysis 

The principal shortcomings of the F-14A bare-airframe flight char¬ 
acteristics at high AOA are wing rock (unstable dutch roll) and roll 
reversal. Otherwise the airplane is quite docile. It is departure- 
resistant, but can be departed if lateral control deflection 1 b sus¬ 
tained without large accompanying rudder. Thus, the following discus¬ 
sion will be devoted to identifying potential factors underlying these 
key differences between the F-14A and F-4J. 

a. Open-Loop Dynamic CharacteriBticn 

A survey plot of root migration with angle of attack for the range 
15 < a < 30 deg and 8 * 0 is shown in Figure 35. This indicates excel¬ 
lent lateral-directional agreement with the previously reported high AOA 
flight characteristics and with the flight traces of Figures 30, 31, and 
33. The dutch roll mode becomes divergent at just under 15 deg AOA and 
stays unstable until approximately 22-23 deg AOA. This is the region of 
reported wing rock. Above 23 deg AOA the dutch roll is stable and well 
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Figure 34, F-14A Model Simulation of 

Flight 236-1 Pullup 












Figure 35. F-14A Root Migrations with AOA 


damped but the coupled roll-subsidence and spiral becomes oscillatory. 
The AOA at which Ng^^ ■ 0 is not shown on this plot because it remains 
positive (as shown in Figure 36). 

The longitudinal short-period mode exhibits relatively low frequency 
and high damping at 15 deg AOA. Frequency then increases and damping 
decreases with increasing AOA until the dutch roll and short-period are 
about equal at 30 deg AOA. Note that this is not in agreement with the 
flight traces for 15 deg AOA (Figure 30), which showed a short-period 
frequency twice that of the dutch roll. 

Referring back to Figure 25, 15 deg AOA is where Cup is a maximum. 
Thus, cross-coupling through Mg may again create pitch frequency cou¬ 
pling and contribute to the nonlinear dutch roll damping, as was the 
case with the F-4J. However for 3^0 the shift in the Figure 35 root 
locations is so small as to be negligible, and this is not the same as 
with the F-4J. 






Ib~ ~20 25 30 

a(deg) 


Figure 36. N 6^y rt Variation with AOA; F-14A 


Eigenvectors for the coupled dutch roll and $hort-*period modes 
of the 6 DOF airframe at « Q ■ 20 deg and 6 0 ■ 1.5 deg are shown in 
Figure 37. There is essentially no longitudinal coupling into the dutch 
roll mode and only alight lateral coupling (mainly rolling motion) into 
the short-period mode. These eigenvectors are noticeably different from 
those for the F-4J (see Figure 14) at a similar flight condition. 

Vector polygons for the dutch roll mode at this flight condition 
are shown in Figure 36. The p equation indicates that £. g and 3f' com¬ 
pletely dominate rolling motion. All damping derives from £ there 
is no cross-coupling influence from All terms in the i equation 

have been doubled to make this polygon legible. The relative size 
of the p end r polygons Indicates that in this AOA region static and 
dynamic yawing force and moment terms are negligible in the dutch roll 
mode. However, the yawing moment that does exist derives primarily from 
Np. The 4 equation has significant magnitude compared to the £ equa¬ 
tion. The key is the dominance of the MglBl term which, as with the 
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4 Equation 



Figure 38, F-14A Vector Polygons for Dutch Roll Mode; 

a 0 “ 20 deg, 2 0 - 1.5 deg 










F-4J, indicates that coupling from the lateral into the longitudinal 
produces the frequency doubling effect. 

Based upon the above argument, it iB concluded that the F-14A wing 
rock 1 b due to a divergent dutch roll which derives from very low 
damping in roll, £ p, rather than lateral-longitudinal cross-coupling 
(since the open-loop lateral-directional dynamic characteristics are 
relatively insensitive to sideslip). The principal static aerodynamic 
cross-coupling is a transfer of energy from rolling-sideslippiug motion 
into pitch motion and this, as with the F-4J, results in frequency 
doubling into pitch and may contribute to amplitude bounding of the wing 
rock. 


Closed-Loop Dynamic Characteristics 


Lateral stick control of roll . Figure 39 presents the key 
parameters for the numerator Nj^. The lower plot indicates that the 
movement of one zero into the RHP with increasing AOA is not as rapid 
for this airframe as for the F-4J (see Figure 17). However, departure 
can occur where one of the real roots is large negative. This plot also 
demonstrates that the roll numerator zeros are quite insensitive to 
sideslip when compared to the F-4J. 

Rudder control of roll . Same characteristics as for the F-4J. 


Stabilator control of pitch . The potential for divergence lies 
6 

with the RHP zero, I/Tq^, in Figure 40 is a plot of l/Tjj^ at 
8 0 « 4 deg over a range of AOA. Note that the parameter scale factor is 
an order of magnitude less than the one used for the F-4J in Figure 21, 
and the largest value (approximately -0.01) is two orders of magnitude 
less than for the F-4J (equal to -1.0). Thus, pilot control of pitch 
attitude should have no noticeable coupling into the lateral-directional 
modes for the F-14A. 
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Figure 40. ^S 8ta ij First-Order Root Migration with ADA, F-14A 


Coupling numerator . The lead coefficient (gain term) for the F-14A 
coupling numerator, ^fi ata ktp» * 8 » ^ or practical purposes, zero. 
Thus closure of either loop (pitch or roll) has essentially no influence 
on the other. 

Summary . The F-14A exhibits somewhat better open-loop static and 
dynamic stability than does the F-4J. It also has considerably less 
lateral-longitudinal static cross-coupling in asymmetric flight, 
and hence little change or degradation in dynamic parameters with 
sideslip. The static aerodynamic cross-coupling derivative values at 
8 0 ■ 4 deg and 15 < o < 30 deg are presented in Figure 41. Comparison 



Figure 41. F-14A Static Aerodynamic Cross-Coupling 

Derivatives 
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of these values with those for the F-4J (Figure 22) shows a drastic 
difference between the two aircraft* 

The most significant maneuver-limiting factor for the F-14A appears 
to be roll reversal associated with movement of one root of the roll 
numerator, into the RHP as AOA is increased. The F-14A exhibits a 

more gradual shift of this root with AOA than does the F-4J. 
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SECTION III 


AEHODTNAHIC CROSS-CODFLING INFLUENCES 


The previous section showed a significant migration of denominator 
and numerator roots for the F-4J due to aerodynamic cross-coupling at 
3 0 - 5.5 deg. In this section literal analytic expressions are devel¬ 
oped to examine in greater detail the influences of specific cross¬ 
coupling coefficients on the airframe characteristic (open-loop) modes 
and the key piloted control numerators for roll due to lateral stick and 
pitch due to stabilator deflections. 

With lateral-longitudinal cross-coupling, at least 5 DOF equations 
of motion (e.g., Figure 13) must be used for analysis. As a result, the 
wealth of existing knowledge on the dynamics of vehicles in three (or 
fewer) degrees of freedom is not adequate, and qualitative insight into 
the motions of the aircraft is sacrificed. 

The important cross-coupling terms in the equations of motion iden¬ 
tified empirically by Porter and Loomis (Reference 29) and Johnston and 
Hogge (Reference 5) are both aerodynamic and kinematic. Effects of 
cross-coupling on the characteristic aircraft inodes were evaluated by 
Hamel (Reference 16) by applying servo analysis techniques to a 5 DOF 
aircraft mathematical model. 

The approach applied next will follow and expand upon that utilised 
by Hamel. 

A. CHARACTERISTICS NODES 

Analysis of an uncoupled (3 = 0) airframe is relatively straight¬ 
forward. Approximations to the classical 3 DOF lateral or longitudinal 
expressions can be developed for a wide range of airframe and flight 
conditiona (e.g., Reference 20). 

However, sideslip characteristically Introduces cross-coupling, both 
kinematic and aerodynamic. Then the only clear simplification one can 
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apply Is the assumption of constant forward velocity, thereby removing 
the velocity equation and associated terms (see Figure 13). This can 
be justified since forward velocity changes comparatively slowly. No 
significant coupling occurs due to velocity, and its effect is primarily 
on the low-frequency characteristics (i.e., the phugold mode). The 
three kinematic equations of Figure 13 can be resolved and substituted 
in the Z, M, and Y equations to reduce the number of equations to five. 

Hamel takes five-degree-of-freedom equations of motion "using body 
axes which have the x-axis initially aligned with the flight path"; 
from the context we interpret that for initial sideslip, the x-axls la 
aligned with the projection of the flight path vector onto the plane of 
symmetry. In expanding the characteristic determinant he separates out 
the effects of the aerodynamic coupling terms N„ and Mg (neglect¬ 
ing any other coupling between the longitudinal set and the lateral- 
directional set) thus: 

A ■ A long A lat + K <M A coupl 
The uncoupled (0 O = 0) denominator terms are 

^long * fl2 + 2C S pUi g pS + m|p 

Alat " (s + 1 /T b )(i + 1/Tr)(s 2 + 2; d u> d s + uj) 

Hamel defines the classical longitudinal (short period) and lateral- 
directional (spiral, roll subsidence, dutch roll) by approximate factors 
(Reference 20). The coupling component, K(0 o )A cou ^, can be manipulated 
to represent a sum of two terms: 

K(Bo>A COU pl - ^B 0 s(s + 1/T 0 )(s 2 + 2c 0 u> 0 s + u> 2 ) 

+ M 0 N'B o (s + 1/TiXs + 1/T 2 ) 

The approximate factors for these coupling terms are reproduced (in the 
notation of this report) in Table 3. 
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TABLE 3. COUPLING NUMERATOR FACTORS 
1/T 0 - -M q 

2 " “^3 “ N r *^o) ®^r 

4 - Ne + YgNj. - (N5/XS)(t£ + *3*r) 

1/Ti - -£p + ( JC.„/N')Cn' - g/U 0 ) 

i/x 2 - ti8/[u 0 C<S/n5)(n; - £;)] 

The roots of the coupling component are evaluated by applying Bervo 
analysis techniques for obtaining the roots of a characteristic equa¬ 
tion, setting it into the form 

1 + M g( N a/*aX s + 1/T l)(» + m 0 

s(s + ]./T 0 )[s* + 2c p u 0 s + to§] 

Figure 42 is a root locus plot for this equation# In this instance 
i Q > 1 and < 0. The "cloped-loop" roots for this case are 

Wo^coupl " <£a«M 8 + *404)(s + .47) 

X [s 2 + 2(-.48)(1.87)s + (1.87) 2 J 

This same technique can be applied on a larger scale to find the 
roots for the coupled charactf '.Stic modes, i.e., the roots of the rela¬ 
tionship 
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Figure 42. Migration of Roots of K(B 0 )A C0U pi with 
Increasing Mg (<* 0 - 21 deg, - 5.5 deg) 


Figure 43 is a root locus plot for this expression showing the migration 
of ta 8p and io d with increasing £ £6 0 . Note that 1/T 8 and 1/T R are 
essentially cancelled by zeros obtained from the approximate factors and 
the preceding "closure." 

The actual 6 DOF roots of the F-4J a Q ■ 21 deg, $ 0 ■ 5.5 deg case 
are shown by the dashed lines and & symbols in Figure 43. There is 
fair agreement between the simplified and exact loci for the dutch roll 
but not for the other modes. The problem, it appears, lies with the 
location of the first-order coupling eeroB. If these were to lie closer 
to the origin, the ui sp and Wj loci would be rotated counterclockwise and 
l/T a and 1/T R would be driven to the right. It is not surprising that 
the approximate factors employed here do not precisely predict the F-4J 
root shift since they were developed in Reference 16 for an airframe 
configuration considerably different from the F-4J. However, they do 
demonstrate the trends and show specifically that: 
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A Actual 6 OOF Roots 



Figure 43. Migration of Roots of Characteristic 
Modes as Cross-Coupling Increases; 
a Q - 21 deg, - 5.5 deg 


• The location of the roots of K (B 0 ^coupl (Fig¬ 
ure 42) and hence the direct Io n of the root shift 
in Figure 43, is controlled by” the sign and mag¬ 
nitude of Mg(N£/£'). 

• T ^ e magnitude of root shift in Figure 43 is con¬ 
trolled by the sign and magnitude of «C. £ 3 0 . 

B. NUMERATOR 

The first step in identifying the Influences of cross-coupling terms 
on the key numerator for lateral piloting control, N$ , require elimi¬ 
nation of all insignificant terms. This is accomplished by a series 
of simplifying steps. The 6 DOF characteristic poles and the zeros 
for the numerator Nj^ are shown in Figure 44, again for the F-4J at 
« 0 ■ 21 deg, 8 q ■ 5.5 deg. Figure 44a represents a case for which the 
cross-coupling terms are zero. At this high AOA the uncoupled numerator 
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Figure 44, N$ stlc /4 Pole-Zero Locations, Coupled vs. Uncoupled, 

F-4J, o 0 - 21 deg, - 5.5 deg 

consists of two real roots, 1/T^ and l/T^* non-minimum phase and 
minimum phase, respectively; and two complex pairs, and u>| p , which 
exactly cancel the uncoupled longitudinal poles. With croBs-coupling 
(Figure 44b) the first-order zeros couple into a complex pair located 
near the divergent dutch roll mode. In addition, the numerator counter¬ 
part to m ap , tUgp, has moved apart from the short-period poles. The 
effects of coupling on the phugoid numerator term, <*>p, are negllble (see 
Figure 44b). This indicates that the velocity equation of Figure 13 may 
be eliminated. The phugoid then reduces to a first-order pole-zero pair 
at the origin, but all other roots are affected only slightly (Figure 
44c). However, even with the terms removed from Figure 13, and 
recognizing that the equations for the Euler angles ^ and 9 contribute 
no significant kinematic coupling and can be deleted, a rather unwieldy 
six-by-six matrix still remains. 

The matrix can be further reduced by use of time vectors evaluated 
for the particular mode of interest. Inspection of Figure 44b shows 
that the ^ zeros of greatest concern are the complex pair which lie 
in the right half-plane, near the already divergent dutch roll poles. 
Therefore further analysis of n| will concentrate on this complex 
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pair, and on the effects of cross-coupling in driving them off the real 
axis, where Figure 44a places the corresponding zeros of the uncoupled 

response. 

For denominator roots the relative magnitude and phase of each 
motion variable in the mode corresponding to any particular real pole 
or pair of complex poles can be determined. The result can be expressed 
numerically as an eigenvector or graphically as phasors. In order to 
determine such characteristics for transfer function numerator roots, a 
pole is driven into each zero to be examined. The poles are driven Into 
the zeros by closing a tight feedback loop, as shown in Figure 45. In 
Figure 44c the real poles are driven into the right half-plane pair of 
complex zeros. Actually the (p/s)^^ transfer function is utilized 
here rather than <|>/6 stlc since the two are almost identical (compare 
Figures 44c and 45). 

Figure 46 presents the resulting phasors for the closed-loop roll- 
spiral mode, oigg, which approximates the numerator "mode" Using 

these phasors, vector polygons have been plotted for this mode (see 
Figure 46) from the equations of motion, Figure 13. Minor vectors 



Figure 45. Closure of p/s ♦ S 8t k t0 Obtain Time Vectors, 
F-4J, a Q - 21 deg, 0 O - 5.5 deg 
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Figure 46. Vector Polygons for Closed-Loop <4* (Approximating c^), F-4J; a D 21 deg, B 0 5.5 deg 



correspond to negligible terms; small polygons indicate little contri¬ 
bution to the modal response. Examination of these vector polygons 
leads to elimination of a number of terms from the equations of Fig¬ 
ure 13. Specifically, the kinematic equations are no longer needed, and 
so the matrix for approximation of the coupled numerator reduces 

to the form shown in Figure 47. The derivatives Y y , and Z w are 

also minimal for this particular airframe and flight condition. How¬ 
ever, under other circumstances they can be quite large, so they have 
been retained to obtain a more general expression for • 

The determinant of this matrix may be developed as a combination of 
uncoupled longitudinal and lateral-directional terms plus cross-coupling 
terms, as follows: 

N ttk : *«atk ' (Z„ + M q + Ma)s - (Ma+ZwM q )] 



This is of the form 

N$ , i [Nx , ) + K . + Km 

stk stk coupled ^ 

where (^^stk) uncoa pi e d consists of two second-order terms: a longitu¬ 
dinal term equivalent to the conventional short-period (denominator) 
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Figure 47. Matrix for Approximation of Coupled Numerator N$ 

approximation; and a lateral-directional term which is the 3 DOF N$ ^ 
(with = 0) (Reference 20). The coupling "gain" terms are: 


and 


■L» f - N ^stk „ 

T* - C0B a o 


stk 


Kn - MgNJ coe «o 


Servo analysis techniques can now be utilized to examine the influ¬ 
ence of each of the cross-coupling derivatives by writing the numerator 
in characteristic-equation form, i.e., 


4 


stk 


^ ( 4 atk ) . . + + % - 0 

stx uncoupled "■ 


and examining the migration of the roots as each of the "gains" is 
varied. 


Setting the coupling term K N to zero, the characteristic equation 
may be written as: 


1 + 




stk uncoupled 


0 
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Figure 48 is a system survey plot showing the effect of varying 

r $ 

positive JL a on the uncoupled numerator The upper plot is a root 
locus and the lower plot a Bode-siggie (see Section II) shoving the 
influence of increasing (negative) gain. These indicate that the com¬ 
plex zero associated with the longitudinal short period, is driven 
to higher frequencies parallel to the jts-axls; and the conventional roll 
zeros, 1/T<j>i and l/T^, are driven apart. The influence of £ £ is found 
to be negligible in this case, but this is due to a very small value of 
the ratio N 6 et j c // £6^ ( at the flight condition being evaluated) which 
is a factor in the "gain" term. 

The prime contributor to cross-coupling influence on the F-4J roll 
numerator is % ■ MgN^ cos a 0 as shown in Figure 49, which is a survey 
for the equation 



where primes have been used to denote that has already been ac¬ 

counted for, and double primes indicate that both K£ and K N have been 
accounted for. K N produces a locus opposite in direction to that uf 
Kjf , since both are pure gains but have opposite sign. The root loca¬ 
tion 1 b much more sensitive to changes In Kjj, as the relative Bode 
amplitudes in Figures 48 and 49 show. 

The "gains" Kjf, and K N are functions of N ti* M 8» N 6gtk’ and 

*^stk* •Co an< * No have keen stressed in the foregoing, it; should 

be obvious that the influence of changes in magnitude of any of these 
derivatives can be equally represented by Figures 48 and 49. Changes 
in sign of any of the derivatives would reverse the direction of root 
migration in either case. 

The results shown here represent the usual derivative signs that 
prevail as AOA increases toward stall, l.e., <£ ^ positive, N' negative, 
and Mp negative. However, Mg depends upon several configuration details 
and can be either positive or negative (Reference 30). For the above 
signs, uncoordinated (8 j 1 0) flight resultB in reduced dutch roll 
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Figure 49. System Survey for Effect of N™ on Ng 
- 3.4251 1/sec 2 ), F-4J; a 0 - 21 deg, g Q = 












stability (Figure 43 or 44) and a corresponding shift of the roll 

numerator zeros. These shifts portend increasing closed-loop roll 

control problems with sideslip (Figure 45). 

The accuracy of the 5 DOF approximation is evident from 

Table 4: the coupled, 6 DOF numerator has been matched almost identi¬ 
cally. With this validation, of the approximation, the following simple 
coupled approximate factors were developed for this flight condition (it 
must be stressed that these literal approximate factors apply to this 

aircraft and flight condition only,and may not be accurate for other 
situations). For the general form 


stk 


4 
2 w 


2C| 


P^sp 
2 


cp 


(s 2 + 25^8 + wjKs 2 + 2;* p u* p s + t«4p) 

-<N£ + Y v ) - /N£M g cos O 0 

Ng cos a 0 + /N£Mg cos a 0 

-(M q + M^) + /NaMp cos a Q 

-M a + Z^q + /N^Mg cos a 0 


Table 4 includes the values given by these approximate factors, 
indicating that they also are accurate and therefore, encompass the key 
derivatives. 

These literal factors clearly indicate that longitudinal-lateral 
cross-coupling in ^ is primarily through the product N^Mg. If 
either derivative is zero, the crosB-coupling terms disappear and the 
approximate factors become similar to the literal expressions for un¬ 
coupled lateral dut-h roll and longitudinal short-period (Reference 20). 
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TABLE 4 


N$ 6tlt NUMERATOR APPROXIMATIONS 


NUMERATOR 



— 

m” 

«P 

6 DOF coupled 

-.531 

.909 

.778 

1.361 

5 DOF (approximate) 
coupled 

-.533 

.892 

.768 

1.370 

Based on literal 
approximate factors 

r- 

.940 

.667 

1.544 


In addition, for non-zero N£ and Mg the aircraft coupling is aensitive 
to sideslip since 


i 3 2 N* . 9 sb 
WZ P ° ~T Z "Set 


so that Figure 49 may be viewed as a locus of root migration with side¬ 
slip for constant Cng a * 

From the standpoint of vehicle design, any airframe changes which 
modify ilther of these cross-coupling derivatives can be expected to 
have a significant influence on high AOA flying qualities and departure 
characteristics. 

C. N^ tab NUMERATOR 

Analysts of the influence of aerodynamic croBS-coupling on the 

0 

numerator for longitudinal control, N 6 sta ^» uses an approach identical 
to that taken with . The first step is to make the development 

tractable by eliminating insignificant terms in the nine-by-nine matrix. 

It wao shown in Subsection ll.A.3.b (Figure 19) that constraining 
velocity does little to change the zeros influenced by cross-coupling, 
i.e., l/Tg^, ^ t R> and utg. Additionally, we find from Figure 50 that the 
equations for the Euler angles 4* and 6 can again be eliminated. This 
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ROLLING MOMENT EQUATION 





Closure, A.q-^ Sstab 

Mode: -f 4 '(~Y g ) p 7 '' - 327 “• 398) 


Figure 50. Time Vectors for Closed-Loop Mode 1/T 
(Approximating 1 /T 03 ), F-4J; 
a = 21 deg, B - 5,5 deg 





still leaves a six-by-six matrix which Is too large and complex to work 
with comfortably. 

Applications of time vectors, as was done with N l at j c » accom¬ 
plished by closing a pitch rate loop with a feedback element of K/s 
(to obtain body-axis 3). Modal response coefficients are then calcu¬ 
lated for the closed-loop pole driven into 1 /Tq g, and the time vectors 
for each equation constructed as shown in Figure 50. Except f^r the 
at equation, few terms are small enough to be eliminated by inspection. 
Further,it seems unwise to restrict the general applicability of the 
results by eliminating the a equation. The main benefit of the time 
vectors is to demonstrate that the modal response for a root at l/Tg^ 
contains considerably more lateral than longitudinal motion. 

Thus, it is necessary to expand the complete six-by-six matrix of 
Figure 51 in literal form, evaluate each of the terms, and eliminate 
products of small numbers. The resulting polynomial is then separated 
into uncoupled and coupled terms to form the approximation: 


N 5 B Ub a «« a tal| { ' ' B W>[* 4 ■ + + + |iCp(Nr + V v )+NiY v + NB Jyn | k 2 

- |of|,Np cos Oo + tfpNfYy + Xiidyn + ( ^r N A “ Bln n a j" 

+ OOB n o + Bln °o|| f 0'p B | s2 " (N r + V» + 7 0 l\] 

+ Niz r s|* 2 + (t~- - Y V - Xpy - piAy + Yv (*P ‘ tan'ajj 


where in the "dynamic axis" system of Reference 18, 

N ^dyn 2 cos “o - *3 8in a o 


®dyn 


cos u 0 


+ Ng sin a c 
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0 

Figure 51. Matrix for Approximation of Coupled Nj 

(♦ 0 - r 0 - 0, 0 - a, |3 0 Small) 8tab 


The terms not containing cross-coupling derivatives are seen to 
separate into a firat-order and quartic function of s. The first- 
order is the conventional literal approximation for the longitudinal 
zero 1/Tg 2> i.e., -Z w . The quartic contains all lateral-directional 
derivatives and, with the exception of the sin a Q terms, is the same as 
the conventional body-axis uncoupled lateral characteristic equation 
(Reference 20), This polynomial thus factors into the parameters 
1/T R , and 1/T S . 

The first cross-coupling term involves as the "gain," a free s, 
and a quadratic containing only yawing derivatives. The second cross- 
coupling expression has N£ as a "gain," a free s, and a quadratic con¬ 
taining only rolling derivatives. 

The kinematic terms Z p and Z r were negligible for N$ but here 
dictate the Influence of the aerodynamic cross-coupling. Ab a result, 
coupling in the 0 * S 8ta ^ numerator is an explicit function of sideslip. 

It Is interesting to note that Mg does not appear in the approximation; 
a quadratic coupling term involving M 6 z 6 B t u b^T 0 was found to be negli- 4 

gible. 1 
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The Influence of the individual cross-coupling terms can again be 
observed by servo analysis methods, setting 


(N§ J - (Nr ) + (N§ J + 

stab coupled stab uncoupled sta ^oda 



is shown in Figure 52. The "poles" 1/T® and 1 /Tq^ couple to create a 
second-order closed-loop root, Wq s , which is on the jui-axis. The dutch 
roll "pole" is moved to the left half-plane, while 1/Tj( is moved slight¬ 
ly to the left. 

Figure 53 shows the additional movement when the cross-coupling due 
to is included. The latter moves things further in the Bame direc¬ 
tion. However, it can be determined that the first-order coupled zero 

in the left half-plane comes from the roll subsidence mode and the zero 

farthest into the right half-plane (previously called I/Tq^) comes from 
a coupling of the spiral mode and the. conventionale 1 /Tq 2 . Since the 
time vectors indicated this mode is dominated by lateral-directional 
motion it is now labeled (l/Tg) . The prime notation in Figures 52 
and 53 1 b used to keep track of the number of "loops" being closed 

and the influence of each "closure." From here on the prime notation 
will be dropped and the lateral-longitudinal coupled zeros identified 
as w®, 1 / t r» etc. 

Since each of the cross-coupling terms has similar influence on 

the migration of zeros, no insight la lost by combining the two coupling 
terras via the concept of the previously noted "dynamic" terra Then 
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The system survey for the combined coupling influence on N fi atab ia shown 
in Figure 54; the coupled roots, of course, are the same. 

Finally, it should be noted that aerodynamic cross-coupling in the 
numerator ^6 sta jj a direct function of sideslip squared . Since JC a is 
proportional to C^^(e 0 ) and N a is proportional to C a p a (3 0 )» the "gain" 
term is 


*«dyn 3 ° ~ (%a B o C0S a ° + ^(Ja 3 ® sin a o^ 

~ C c a 3a . c° 8 «o + c ng a 8in a o)^o 

Thus the zero 1 /t| is drastically influenced by static aerodynamic 
cross-coupling and magnitude of sideslip. 

D. SUMMARY 

Literal analytic expressions have been developed to show the influ¬ 
ence of static aerodynamic cross-coupling derivatives and asymmetric 
flight on denominator and numerator roots. The analysis has shown that 
the shift in denominator roots depends upon all three static cross- 
coupling derivatives, NJ, and Mg. The direction of root migration 

is controlled by the sign and magnitude of MgN a / The magnitude of 

root migration is determined by ^ a 8 0 . 
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Figure 54, System Survey for Effect of Combined Cross-Coupling 
on N$ stab ; F-4J; a G -= 21 deg, $ 0 “ 5.5 deg 
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For the F-4J the. principal croBB-coupling Influence on the roll 
numerator, i® due to the term N^Mg cob a Q . The literal approxi¬ 
mate factors which describe the coupled numerator roots are 


■ -<NJ + Y v ) - /N'Mg cos « 0 

2 • + - 
10 $ ■ Ng cos a Q + /cos ot Q 


2cJ p 4p " -<Mq + M$) + /N£Mg cos ct 0 
“sp * " M a + ZwMq + ^ K a M g 008 a o 

The pitch numerator, N 6 8ta jj» Influenced by a more complicated 
combination of kinematic and aerodynamic terms. Existence of a first- 
order R'rfP zero depends upon 

£ q Ng — <C gNg — Yy( if r COS <X 0 ** p Bln a 0 ) 

oCq cos qIq + sin Uq r 

The above approximations may not apply to configurations significantly 
different from the F-4J. 
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SECTION IV 


CONFIGURATIONS SELECTED FOR SIMULATION 


The seven aerodynamic derivatives which have been shown to dominate 
open- and closed-loop dynamic characteristics at high AOA are summarized 
in Table S. These are comprised of six static stability and cross- 
coupling derivatives and one dynamic damping derivative. Together they 
limit the useful flight envelope via roll reversal and wing rock tenden¬ 
cies, and they limit the safe flight envelope via departure suscepti¬ 
bility. One goal of the piloted simulation was to demonstrate this 
premise. 


TABLE 5 


KEY MANEUVER-LIMITING PARAMETERS 




Negative Ns a or N$ D 

Roll reversal 

Positive Mg 

Pitch up 

tb Xp ) 

l Wing rock 

Ng. N' a > 

\ Nose slice 


\ Roll divergence 


This section presents the vehicle aerodynamic configurations se¬ 
lected for the piloted simulation and the predicted high AOA dynamic 
characteristics. 

Four "configurations" were obtained by altering the key aerodynamic 
coefficients of a baseline. The three altered configurations incor¬ 
porated changes in, respectively, C£p(a), C$(a,B), and C n (a,(J) combined 
with a change in C m g, to provide specific alteration of key open- and 
closed-loop parameters discussed in the preceding two sections. The 
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attempt was to Influence the denominator only, the numerator only, and 
the combined denominator and numerator, thereby to provide distinct 
variations in high AOA departure susceptibility, onset warning, and 
motion severity. The influence of artificially altering aerodynamic 
coefficients via a lateral-directional augmentation system was also 
investigated with two of the configurations. 

A. AERODYNAMIC CONFIGURATIONS 


1■ Configuration A 


The baseline case is iden 'led as Configuration A. This F-4J- 
based configuration has been thoroughly reviewed in the preceding sec¬ 
tions. The key open- and cloaed-loop dynamic parameters are summarized 
in Figure 55 for comparison with the other configurations discussed 


subsequently. Dutch roll and numerator root migrations are shown for 
increasing a at zero 0 (solid lines) and at 5.5 deg 8 (dashed lines). 


The solid triangle (A) indicates the AOA above which C n ^ becomes 
negative at zero 8. Recalling the dlBcusBion in Section II about the 


shortcomings of linearized frozen point analysis, the important aspects 
t.o note here are: 



The low dutch roll damping at AOA below C 
(wing rock). 


n 3dyn 


0 


2 

or 1/T^ negative at AOA above 18 deg or at 
small 3 (roll reversal). 


e 1/T*. an< * l/Ts both large and negative when 6 j 4 0 
(nose slice parameters). 


e Dutch roll and roll numerator roots sensitive to 
sideslip. 

The plots in the lower left corner indicate the magnitude of the 
aerodynamic coupling terms evaluated at 5.5 deg 8. In the AOA range 
from 20 to 30 deg, is cf approximately the same magnitude as but 
is of opposite sign. In this same AOA range, is of about the same 
magnitude as Njj and of the same sign. For this relatively mild sideslip 
condition (5.5 deg), the static cross-coupling derivatives are of equal 
magnitude to the basic static derivatives. 
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2. Configuration B 


For Configuration B the roll rate damping coefficient, Cj^, was made 
large, negative and essentially invariant over the range 10 < a < 35 deg 
(see ^sute 56). The object of this change is to increase the 
dutch roll damping and hence reduce the wing rock tendency. Since all 
other parameters remain unchanged, this should not Alter the basic roll 
reversal tendency or the nose slice departure characteristics; but It 
should eliminate or reduce wing rock as a warning of impending depar¬ 
ture. However, note again in Figure 57 the sensitivity of the dutch 
roll mode to sideslip. At 8 - 5.5 deg there is little difference 
between Configurations A (Figure 55) and B (Figure 57), which indicates 
that the Influence of «Cp is small with respect to that of jCjj. 



Figure 56. 


Change in Co 

P 


for Configuration B 
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3. Configuration C 


For Configuration C the roll due to sideslip coefficient, C$g, was 
altered (see Figure 58) to approximate that of the F-14A over the range 
15 < a < 45 deg, taking into account the difference in roll inertias of 
the two aircraft, i.e., 


%2 


7-14 


r *F-4 

1j «F-14 


Figure 59 presents the parameter plots for Configuration C. The 
upper left figure shows that the change Introduced in Cjj,^ and Cg a has 
shifted the dutch roll mode toward the left. It remains in the left 
half-plane for 8*0 and only bulges into the right half-plane over a 
small a range at 8 * 5.5 deg. A major difference in this configuration 
occurs in the roll subsidence and spiral modes, which in the previous 



Figure 58. Change in for Configuration C 
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cases have remained on the real axis and therefore have not been shown. 
In this case, 1/T R and i/T 8 couple to form a "lateral phugoid" mode, 
u f)R , as AOA increases. At 36 deg AOA the frequency of this mode is 
approximately twice that of the dutch roll, with approximately the same 
damping. Comparison of Figures 35, 39, and 59 show this to be a valid 
representation of F-14A root locations, but mode identity may be differ¬ 
ent. 


Comparison indicates that increasing 6 shifts the uig R locus toward 
the right half-plane. However, both the dutch roll and lateral phugoid 
modes are now relatively unaffected by 3 . 

The upper right-hand plot for is essentially the same for Con¬ 
figuration C as that for Configurations A and B. While this configura¬ 
tion has better open-loop stability, the closed-loop characteristics 
should be little different from the previous two configurations because 
the roll numerator roots are essentially unchanged. 


Absence of the solid triangle symbol indicates that 
tive throughout the AOA range. 


Sdyn 


is posi- 


The lower left-hand figure shows that the term (which derives 
from the variation of Cj^ with «) is small and negative. Since this 
configuration also has a large increase in £3 (see Figure 58), the 
aerodynamic cross-coupling derivative is now much smaller than £ | 
throughout the angle-of-attack range. Comparing the values of £ £ for 
Configurations A and C with the differences in denominator root migra¬ 
tions with 8 makes it apparent that /j is the dominant purameter in 
shifting the Configuration A dutch roll mode into the right half-plane 
when 8^0. The lower right-hand plot of Figure 59 indicates that the 
change in and £ j$ has decreased 1 /t| at AOA above 20 deg. 
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4. Configuration D 


For the fourth case the yaw due to sideslip coefficient, C n g, was 
modified to approximate that of the F-14 as shown by the dashed line In 
Figure 60, In addition, for this case only, the sign of C m g was re- 
versed to provide positive pitch with sideslip. As predicted by the 
approximate factors of Section III, the result Is an airframe quite 
Insensitive to sideslip (Figure 61) and, most pertinently, having numer¬ 
ator roots (l/T^ and 1/T^) which do not penetrate as far into the right 
half-plane. 

Roll reversal is delayed until a > 20 deg and 
tive until a * 25 deg. 

aaaao 
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Figure 60. Change in C n ^ for Configuration D 



C n ^ remains posi- 








5. Anticipated High AQA Maneuver- 

Limiting Characteristics 

Review of the preceding four configurations should go a long way 
toward sorting out the relative Importance of denominator (open-loop) 
vs* numerator (closed-loop) characteristics on high AOA departure sus¬ 
ceptibility and severity. Configurations A, B, and C have nearly iden¬ 
tical RHP numerator root locations and therefore should have similar 
closed-loop divergence tendencies, although the motion (i.e., roll vs. 
yaw) may differ somewhat due to Cjj, pl Cg^, and Cj, a changes. On the other 
hand, the Configuration D numerator roots do not penetrate so far into 
the RHP and are relatively insensitive to sideallp. Thus one might 
expect this vehicle to appear less susceptible to or have less severe 
departure. 

Alternately, if departure susceptibility or tendency should turn out 
to be strongly influenced by dutch roll (open-loop) root location, then 
one would expect Configurations C and D to be similar Bince both are 
near neutral stability and relatively insensitive to sideslip. Any 
differences between Configurations A and B would be dependent upon side¬ 
slip. 

The dynamic stability parameter, C na , has been referred to pre- 

"dyn 

vlously. The AOA at which this parameter becomes negative was identi¬ 
fied in Figures 55, 57, and 61. Figure 62 presents plots of 

vs. a for the four configurations. Based upon the criterion 
> 0.002 for departure resistance (Reference 31), then only 

, 1 

Configuration C should be resistant, Configurations A and B should be 
quite departure prone, and Configuration D might be mildly departure 
prone. 

Another criterion, Reference 32, employs both open- and closed-loop 
parameters, i.e., 


n ^dyn 

of 
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Figure 62. C na for Configurations Simulated 
s dyn 


For a number of years it has been recognized (Reference 33) that unde¬ 
sirable (w^/wjj) ratios lead to roll reversal and pilot-aggravated wing 
rock (PIO). However, in Reference 32 Weissman plots C n& vfl * LCDP and 
identifies four regions of increasing departure and spin Buscepticility 
and severity (Figure 63) based upon experience with a number of airframe 
configurations. The loci of our unaugmented airframe parameter values 
over the range of 16 < a < 35 deg are shown in Figure 63; circles repre¬ 
sent Configurations A and B (since C* p has no influence on either param¬ 
eter), squares represent Configuration C, and triangles Configuration 
D. The criteria predict high departure/spin susceptibility with strong 
rolling departures for Configurations A and B, low spin susceptibility 
and mild rolling departures for Configuration C, and no departures for 
Configuration D. Thus, all regions of departure susceptiblity and 
severity are exercised with the aerodynamics selected. 

The Figure 63 prediction that Configuration D should be less suscep¬ 
tible to departure is in agreement with the first assessment above based 
on numerator root locations. 


: 


1 
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Figure 63. Weisstaan Criterion Predictions for 
the Simulation Configurations 








B. AUGMENTED CONFIGURATIONS 


Lateral-directional augmentation was synthesized for Configuration A 
(the most departure prone) and for Configuration C (C n ^ always posi¬ 
tive). The augmented configurations are identified as A 2 and C 2 , re¬ 
spectively, while the unaugmented configurations will be referred to as 
A^ and Cj. The augmentation mechanization is representative of that 
used in current high-performance fighters to improve both open- and 
closed-loop handling parameters, i.e., washed-out yaw rate and lateral 
acceleration feedbacks to rudder to increase dutch roll damping; roll 
rate feedback to differential stabilator to reduce or eliminate wing 
rock and augment dutch roll damping; and lateral stick to rudder cross¬ 
feed (SRI) to reduce or eliminate adverse aileron yaw for pilot inputs, 
thereby providing a positive shift in the closed-loop parameters LCDP, 

«|» and 


Undesirably, the roll damper also produces adverse yaw when opposing 
uncommanded roll; it may increase nose slice departure tendencies and 
pro-spin surface deflections. A block diagram for the augmentation 
mechanization is shown in Section V, 

The root migration surveys for Configuration A 2 are shown in Fig¬ 
ure 64, Comparison with Figure 55 shows that at 0 - 0 the dutch roll 
damping is increased considerably; it does not go unstable until AOA 

exceeds that for C n „ becoming negative (A symbol). Also, the AOA at 
? B dyn 

which becomes negative has been increased approximately 5 dcg and the 

negative values of both l/T^ and 1/T® have been reduced. Thus, within 
the augmentation system authority limits and near zero 0, Configura¬ 
tion A 2 stability is increased and departure susceptibility decreased up 
to about 23 deg AOA. However, for 0 j 4 0 the augmentation benefits are 
overcome and the dynamic parameters revert back to air 4 rarae-alone 
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Figure 64. Key Parameters for Configuration A 2 


Root migration survey plots for Configuration C 2 are presented in 
Figure 65. Comparison with Figure 57 again shows Improvement in dutch 
roll stability. However, the most dramatic change is in ui| which now 
remains positive up to 35 deg AOA (but note is small and negative 

beyond 31 deg AOA). Thus roll reversal problems are essentially elimin¬ 
ated and numerator root sensitivity to sideslip is decreased consider¬ 
ably — again subject to augmentation system authority limits. 

The change in departure and spin characteristics cannot be pre¬ 
dicted for these configurations using the static aerodynamic parameters 
of Figure 63 because the augmentation system influence is frequency 
dependent. But a somewhat similar plot can be constructed using the 
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Figure 654 Key Parameters for Configuration C 2 


2 * 

parameters and C^d* Such a plot is presented in Figure 66 for 

both unaugmented and augmented versions of Configurations A and Ct The 
2 

Intercept of the axis with has been arbitrarily shifted to en¬ 

hance similarity for the unaugmented configurations with Figure 63, 
Note that the root plots for the unaugmented Configurations Aj and 
(solid lines) are quite similar to those of Figure 63. The plots 

(Figure 66) for Configurations A 2 and C 2 (dashed lines) are moved fur¬ 
ther into the upper right quadrant, which is identified as departure 
resistant in Figure 63. Similar interpretation would forecast A 2 to 
be departure resistant up to about 24 deg AOA but to have moderate to 
severe departure tendencies at higher AOA, while C 2 should have no 
departure tendencies. 


Despite the direct relationship between C n and it is move 


appropriate to plot C^u)^ because the chauges in 
of C n and are in closer agreement. 

8 dyn 


dyn 

sigr 


gn (or lack thereof) 


105 











2 

Figure 66. Versus for Augmented and Unaugmented 

Configurations A and C 


C. SUMMARY OF CONFIGURATIONS AND 

PREDICTED HIGH a CHARACTERISTICS 

Table 6 summarizes the six-configuration matrix employed in the 
piloted simulation. For A^ a sequence of roll reversal, wing rock, nose 
slice, and finally strong rolling departure with high spin suscepti¬ 
bility is predicted with increasing AOA. For A 2 the augmentation system 
and crossfeed are expected to minimize or eliminate the roll reversal 
and wing rock predeparture warnings. The predicted characteristics are 
nose slice followed by moderate rolling departure. For B, predicted 
characteristics are roll reversal, nose slice, and strong rolling depar¬ 
ture. Warning In the form of wing rock should not be present due to the 
large roll damping of this configuration. Note that Configurations A 2 
and B allow comparison of high AOA stall/departure characteristics with 
high aerodynamic roll damping (Configuration B) and with artificially 
augmented roll damping (Configuration A 2 ). 

The unagumented Configuration C^ is predicted to exhibit roll rever¬ 
sal, wing rock, and mild rolling divergence with increasing AOA but to 
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TABLE 6. CONFIGURATION MATRIX 


CONFIGU¬ 

RATION 

FCS 

AERODYNAMIC VARIANT 

PREDICTED 

CHARACTERISTIC 

A 1 

Basic 

Basic F-4J 

Wing rock 

Roll reversal 

Nose slice 

Roil departure 

A 2 

Augmented 

Nose slice 

Roll departure 

B 

Basic 

Increased Co 

P 

Roll reversal 

Nose slice 

Roll departure 

C 1 

Basic 

Increased Co. ) 

p [ 15 < a < 45 

Decreased C« ) 

*a 

Wing rock 

Roll reversal 

Roll departure 

c 2 

Augmented 

9n«. > 0 

None ’’dyn 

LCDP > 0 

D 

Basic 

Increased C n „ ) 

M o > 15 

Decreased C n ) 

“a 

Positive C 0 g 

Wing rock 

Roll reversal 

Pitch up 


have no nose slice. Augmentation is employed with C 2 to Improve maneu- 
vering control and to determine if it will degrade the departure resis¬ 
tance of this configuration. On the basis of C nfl and (or LCDP), 

Pdyn * 

no departure tendency should be anticipated for C 2 . Finally, D should 
exhibit roll reversal and wing rock warnings but little lateral-direc¬ 
tional departure tendency. However, the change in sign of C,^ to pro¬ 
vide positive pitching moment with sideslip should result in pitch-up, 
which would be expected to aggravate any high AOA departure tendency. 
All other configurations have negative C^, should pitch down with 
increasing sideslip, and hence should require more effort of the pilot 
to maintain high AOA (longitudinal stick cue). 






SECTION V 


PILOTED SMUTATION 


The purposes of the piloted simulation were to validate the influ¬ 
ences of the key aerodynamic coefficients In determining departure 
characteristics, evaluate the Influences of varied maneuver limiting 
factors on high AOA maneuvering control, and identify potential flying 
qualities criteria. While aircraft motion is highly desirable for such 
validation and assessment, departure and post-stall gyration (PSG) 
severity obtained with the 6 DOF analytic models indicated that a 
moving-base simulation would offer very little benefit. Nose slice 
motion was so rapid that 130 dcg/sec yaw rate could be reached within 
two seconds after departure onset; peak angular accelerations as large 
as 

p ■ 500 deg/sec^ 

q ■ 200 deg/sec^ 

i - 150 deg/sec^ 

were obtained. Washout requirements to prevent hitting displacement 
stops would have to be so rapid as to negate motion benefits. There¬ 
fore, a fixed-base simulation was selected having an unrestricted out- 
of-cockplt horizon display capable of such rates and accelerations. The 
aerodynamic model was also expanded to accommodate Inverted flight 
maneuvers. 

This section presents an overview of the simulator, including dis¬ 
plays, airframe, and flight control system modeling, cockpit, and force 
and moment equations. 



A. SIMULATOR 


The simulation was performed in the McDonnell Aircraft Company's 
fixed-base Manned Air Combat Simulator dome, a 20 ft hemisphere identi¬ 
fied as MACS-1. Physical aspects of the simulation are summarized in 
Figure 67. The horizon and target are projected on the inside of the 
hemisphere. The. cockpit is located at the center of the dome. The out- 
the-window, head-up (HUD) and head-down (HDD) displays, and cockpit 
layou: are as indicated in Figure 67. Seat cues consisted of normal 
acceleration and buffet motion provided through an inflatable seat 
bladder. The frequency and amplitude of the buffet oscillation varied 
as the aircraft angle of attack increased, starting at about 14 deg. 

A TV, projection of a gimballed model target aircraft provided a 
maneuvering tracking task. Target motion equations provided realistic 
maneuvers while allowing the aircraft to fly well beyond the departure 
region for the test aircraft. The target was controlled by the computer 
operator via a special control panel. 

B. HEAD-UP DISPLAY 

The head-up display (HUD) provided Information on aircraft attitude, 
heading, airspeed, altitude, and trim conditions. Figure 68 illustrates 
the HUD presentation as originally installed. Display of normal accel¬ 
eration provided a reference for trimming longitudinal control to 1 g 
prior to the start of a run. Trim thruBt was set by adjusting throttles 
until change from trim thrust registered near zero. The Flight Path 
Angle Ladder displayed both flight path elevation angle and roll atti¬ 
tude. Orientation of the total velocity vector was provided through 
the Velocity Vector Indicator; when the velocity vector was outside the 
HUD field of view, the Indicator remained at the edge of the HUD and 
blinked. 

At pilot request, the HUD was simplified late in the simulation 
program. Airspeed, Heading, and Altitude bar displays were removed. 

The fixed reticle tracking gunsight was a part of the HUD. 
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Figure 68. Head-Up Display (HUD) Used In Simulation 


C. FLIGHT OOHTBOL SYSTEM 

The manual flight control mechanization (stick, pedals, surface de¬ 
flections and rates, etc.) represented the basic F-4 aircraft (e.g., see 
Figure 2). 

A special lateral-directional stability and command augmentation 
system (SCAS) was mechanized as shown in Figure 69. The SCAS feedback 
gains are production F-4 values. The stlck-to-rudder crossfeed gain 
and shaping were optimized for 10 < a < 20 deg and faded to zero at 
a < 10 deg to prevent adversely influencing the low AOA handling quali¬ 
ties. 

D. EQUATIONS OF MOTION 

Slightly modified versions of the standard McDonnell Aircraft Com¬ 
pany piloted simulation equations of motion were used. These are essen¬ 
tially the same as those presented in Appendix I, Part III, i.e., moment 
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equations in body centerline axes, force equations in wind axes, and 
aircraft orientation in standard Euler axes. The only change required 
was to expand the small angle a and 3 assumptions and to eliminate 

the tan 0 and cos"' 1 ’ 0 expressions of the <1 equation which gave discon¬ 
tinuity problems as 0 approached 90 deg. This required substituting 

the <requation for 4, i.e., 

4 - q - tan 0(p cos a + r sin a) + Zy/CmVip cos 0) 

in ■ qVj cos 0 - V T sin 0(p cos a + r sin a) + 2^/m 

and using the simple expressions a ■ sin - * (w/V^,) and 0 ■ sin - * (v/V T ) 
for entry.into the aerodynamic look-up tables. 

B. AERODYNAMIC MODELS 

Since the major interest in the simulation was the analysis of 

departure onset, particular attention was given to realistically model¬ 
ing aerodynamic data for the range of 0 < a < 45 deg and 0 < ±30 deg. 
The basic aerodynamic coefficients were selected with the goal of pro¬ 
viding a reasonably accurate dynamic model of the F-4J over this a and 0 
range. To allow greater aircraft motion freedom typical of severe 
departures, post-stall gyrations, and spins without encountering unreal¬ 
istic discontinuities in the data, certain approximations were made for 
a > 45 deg, a < 0 deg, and 0 > ±30 deg, i.e., 

• For a > 45 deg the coefficients were faired to 
110 deg, then extended to 180 deg by either main¬ 
taining the value at 110 deg or fading the coef¬ 
ficient to zero at a • 180 deg. 

• The coefficients were assumed to be similar in 
value for positive and negative A0A so that all 
derivatives were taken to be even functions of A0A 
(i.e., mirror images about a - 0) and the coeffi¬ 
cients and C m to be odd functions of AOA. 

• Variation of the coefficlonts with 0 was assumed 
to be linear for 0 > 30 deg. 

The resulting aerodynamic, coefficients are shown in Figure 70. 
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It is recognized that the assumptions made above could result in 
considerable error for extreme angles of attack and sideslip; however, 
the use of full-ranged data allowed the simulation to be operated 
throughout a violent departure without demanding an immediate return to 
initial conditions when lower aerodynamic data limits were reached# 

The aerodynamic data were stored in the digital computer as look-up 
tables# Total forces and moments were computed digitally and fed to the 
analog systems of the simulator to drive the cockpit displays. 

Physical dimensions of the airplane were based upon the P-4J. Iner¬ 
tias and weights are summarized in Section II# 

F. DATA RECORDED 

Parametric data were recorded via three eight-channel Brush record¬ 
ers. All Euler angles, body axis rates, and control deflections were 
recorded, as well as accelerations, altitude, velocity, and thrust set¬ 
tings. 

For a limited number of runs the total body axis accelerations 
(p» q. r) were recorded and inertial and aerodynamic components of 4 
were also recorded separately. 

G. TASKS AMD MANEUVERS 

Piloting tasks were selected which exercised the open- and closed 
loop departure parameters of Interest. These were divided into two 
phases. The first consisted of familiarization maneuvers typical of 
"feeling out" stall/departure motions, warning, etc., of a new airframe. 
These are modeled after the suggested test matrix of Reference 34 plus 
an additional aggravated input maneuver for determining departure sus¬ 
ceptibility, proposed in Reference 33. The second phase consisted ot 
tracking tasks and maneuvers as suggested in Reference 36. 
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1. Handling Familiarization and Aaaameut 


Five basic familiarization (F) maneuvers were employed: 

• Straight-ahead stall (Fj): 1 g , fixed thrust, 

longitudinal stick pull to produce d * 1 deg/sec 
until stall departure. 

• Constant attitude stall (F?): 1 g, wings level, 

holding constant pitch attitude while slowly re¬ 
ducing thrust until stall/departure. 

• dank-to-bauk turns (F 3 ): Constant altitude, 

±60 deg bank-to-bank turns at increasing AOA until 
departure; rudder coordination optional. 

• Wind-up turn (F^): Constant thrust, wind-up turn 
to stall departure. 

• Full-stick-deflection maneuver (F^)i From 60 deg 
bank, rapidly ramp full aft longitudinal stick 
followed by full lateral stick out of the bank. 

After 8 sec neutralize controls. 

With the exception of maneuver F 5 there was no restriction on the 
type or magnitude of conttole used by the pilots. In fact, pilots were 
advised to investigate various control techniques Just as they would in 
a stall evaluation of a new aircraft. Maneuver Fg (from Reference 35) 
was evaluated as an indicator of departure susceptibility. 

Hm Precision Tracking Evaluation 

The target aircraft (a TV picture of a model airplane, projected 
onto the dome at the proper aspect and size) was programmed to perform 
three basic maneuvers during a slow pull-up. Target climb rate was 
fixed to produce 4*1 deg/sec to guide the subject pilots into repeat¬ 
able tracking runs of approximately 30 sec minimum duration before 
departure. The tasks started with the subject aircraft trimmed at 
u * 12 dig, essentially in trail. The pilot attempted to stay in the 
target's 6 o'clock position as if trying to overextend a firing 
solution. The tracking (T) maneuvers were: 
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• Straight-ahead pullup (T^): Conatant-thrust, 

straight-ahead pullup keeping pipper on target 
until aircraft no longer controllable. 

• Climbing roll-reversals <T 2 >: Track target air¬ 
craft through series of climbing roll rever¬ 
sals ($ * ±40 deg) until aircraft no longer con¬ 
trollable . Constant'thrust. 

• Wind-up turn (Tj): Starting wings-level, track 
target into constant-thrust 60 deg bank wind-up 
turn until aircraft' no longer controllable. 


3. Departure Recovery 

Recovery techniques employed were at the discretion of the pilots, 
based upon their experiences in fighter aircraft. Recoveries were ini¬ 
tiated either: 

• Immediately upon positive indication of uncora- 
manded motion <R^), or 

• After about a 3 aec delay (R 2 )» 

Recoveries were also attempted In the SCAS-on configurations (Cases A 2 
and C 2 ) by turning the SCAS off upon departure (R 3 ). 

0. FLYING QUALITY RATING SCALKS 


Previous stall/departure simulations for flying quality assessment 
(e.g., Reference 5) have shown the familiar Cooper-Harper scale (Fig- 
w-e 71) to be Inadequate because departure ia centered on the first step 
i v the decision tree, i.e., controllability Instead of performance. It 
was preordained for this simulation that control would be lost (CH ■ 10) 
because the configurations were specifically selected to provide differ¬ 
ing departure warning, severity, and recovery In order to identify gra¬ 
dations in these attributes. Therefore, one requirement of this program 
was to develop and validate a morn appropriate stall/departure flying 
quality rating scale. 

The Reference 37 experimental program succeeded In developing rating 
scales for a similar loss-of-control and recovery situation induced 
by flying into a very strong wake vortex. This experience indicated 
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that each attribute must be evaluated individually and then the overall 
hazard assessed. Figure 72 is a tentative departure rating scale adapt¬ 
ed from the wake vortex experience. Here the stall/departure attributes 
are divided into four open-ended scales. The first is departure onset 
warning — clarity or lack of warning. Next, the attitude and attitude 
rate excursions are evaluated. Recovery characteristics are separately 
evaluated (whether the recovery was rapid and achieved by simple control 
application or whether it degraded down to the point of no recovery). 
Finally, an overall assessment of hazard 1 b rated. 

In anticipation that overall hazard assessment (acceptability) might 
vary with aircraft primary mission, provision is made for indicating 
scalar differences between training and operational squadron usage. For 
example, an operational aircraft might have very clear warning of depar¬ 
ture onset and extreme excursions if the pilot persisted, and be diffi¬ 
cult to recover, but be given a relatively good hazard rating, i.e., a 
2, on the basis of the very clear and consistent warning. Similarly, a 
training aircraft might have little warning of departure onset and rela¬ 
tively large attitude excursions, but recover very easily and quickly by 
merely releasing the controls, and again might have a relatively good 
overall hazard rating, such as a 2, on the basis of recovery. 

As an adjunct to the rating scales, the pilot is also asked to pro¬ 
vide qualitative information as to the nature of the aircraft motions 
observed. 

Although the rating scale of Figure 72 was conceived to cover both 
open- and closed-loop situations and specifically is not spin oriented, 
it was employed in a separate departure/spin susceptibility piloted 
simulation (Reference 38) to obtain an early assessment of its useful¬ 
ness. Results were somewhat inconclusive since the Reference 38 task 
involved preprogrammed, open-loop pilot inputs to assess the departure/ 
spin susceptibility of an airframe having variations in surface author¬ 
ity and actuator rate limiting characteristics. Nevertheless, valuable 
comments were obtained (Reference 39) along with a recommended expansion 
and restructuring of the pilots' qualitative assessment to provide 
greater flexibility and provoke additional commentary to identify what 







DEPARTURE ATTITUDE AND ATTITUDE RECOVERY OVERALL 
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the pilot is rating. Accordingly, a loas-of-control/departure/recovery 
debriefing guide was prepared, Table 7, to augment the rating chart. 
This guide addresses the four topics of the Figure 72 rating scale plus 
an additional section devoted to post stall gyrations (PSG) characteris¬ 
tics and severity. Provision is toads under each topic for the numerical 
rating (STI) from the Figure 72 scale. The pilot is requested to answer 
those questions appropriate for the maneuver flown. Most questions per¬ 
tain to characteristics of the vehicle motion and are constructed to 
encourage simple yes/no, multiple choice, or short written answers. 
Cooper-Harper (CH) ratings are also requested for the two tasks which 
Involve pilot effort to exert control on the vehicle (l.e., prevent 
departure and accomplish recovery). Finally, it should be noted that 
the questionnaire was prepared for use in flight as well as simulation. 
Thus, some questions are not appropriate for fixed-base simulation 
(e.g., assessment of vehicle acceleration). 




TABLE 7. LOSS OF CONTROL/DEPARTURE/RECOVERY 
DEBRIEFING GUIDE 1 . 


a) U warning cloor/uncUAr/mmtKiattnif __ 

b) What U tht nalurt of tilt warning? (Aircr.Wt motion*, 

vibration*, inatrlnn.’it iniilcotIcmi, visual cu**, moUirn cau», cant 
iyatcm (oil cut*} , 

e) la warning mankbd by bomb othar aircraft ciiaracturiati 
d) 1 b tht wmning luch that It allow* th# aircraft to bw 

flown elocsr to the Unit of tha anvalop#? ... __ 

i) Do«« cha warning Intiifara with tno pnn.uy tas.i or 

raduca minion of factlvtnan? . __ 

0 Oo th* aircraft motion*/vibration* cunotltuu j mmm 

huacd worthy of 4 flight rootrictionf __ . . _ 

i) ti tin warning to inada^uau that .t ilight ro»t tut ion 
would b« nacuiitiy to prtvonc lot* of control? _ _ __ 


b; Wivii wot* ilia aircraft motion*’/ (Wing rock, no 1 
pitch up, roiling daparsurt, dlvtrgarU one 1 1l,uInna) 

a) What utt th* atvarlty of chv aircraft motion? _ 
Mild - Mild Arbitration two* 

Incurnadiati - Rapid motion or bucaitratluii in ■ 

mote axis* 

Savore- Very rapid notion or uccularatlon In i'm 


1) How large wort th« changon in alrci'jil •iUUuJu?. 

2) How (bat war* tht ratta? _ . _ __ 

J) How largb uara tha tcsalarailona? ___ _ 

d) Ware th* aircraft motion* and attuudaa roLjfjn.i.ibluI. 
a) Wtra thb aircraft motion* diaori*ntlug and/or del'll it 
l) According to Mll-S-flJMi lUfinltiuna, huw would the r 
'.unet/iuicapclbiUty to daparturt b* doacribud for ^ivl* manuuvur? 
fi) Should pilot action praline or daljy do|wimiro? (U n 
pilot action takon io atato) , __ 

1) Wliuc pilot actiona uuta takon?__ 

2) uiiat uara/wauld bb tin* demand* on thb pilot to 

prtvbtu dapartura? (Kafur to Coopar-llarpar Seal*) _ . 

3) Old pilot actiona agyravato dboartura? 


b; hliat wa* thb Bbviriey of tin a (rural t notion? 
Mild - Hild anublar.itlona and ratta 
Intarmartlata - Rapid motion or auualaratiun 1 


Savar# - Vary rapid motion i 
l) W«rt thb dhana** In air 


dl»or Writing .iml/or dibil 


ptrloil of Lima, short 


Ag|;r.iviHl.ij| - to th-a ill-of •••out r 
d) Wh.it wor« the dnniii.il* c>i til* pi ly 
ocovary? (Rof*r cu Coopor-IUrpor Snl* when 

••fda dns 'rllting danaml* mi pLlm, Mmli'v •«. 
unaldurlnu that thii uik began with an mil-■ 
V ivarall llaiarl r.rh*-«r 
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SECTION VI 


SIGH AOA FLYING QUALITIES ASSESSMENTS 

Effects of the maneuver-limiting factors on high-angle-of-attack 
flying qualities were assessed by three pilots utilising the prescribed 
set of non-tracking and tracking maneuvers. The two principal pilots 
were highly qualified U.S. Air Force fighter test pilots, graduates of 
the Air Force Test Pilot School at Edwards AFB, with a broad range of 
practical experience. The first pilot, RC, has flown a varied selection 
of aircraft, from the F-4 to the B-52; the second pilot, JF, was experi¬ 
enced in operational fighter and attack aircraft, including the A-9 
stall/spin evaluation, and is a former instructor at the Test Pilot 
School responsible for stall/spin curricula. The third pilot, RH, was a 
handling quality engineer/general aviation pilot who also accomplished 
all of the Initial simulation debugging, procedural checkout, primary 
pilot brieflngB and debriefings, etc. 

The combination of airframe configurations, familiarization and 
tracking maneuvers, and recoveries resulted in the 98-cell matrix of 
Table 8. This matrix was flown first as a sequence of runs progressing 
horizontally across the various maneuvers for a given airframe. After 
the complete matrix had been accomplished once, a second series was per¬ 
formed progressing vertically through the airframe configurations for 
each tracking tank. 

A total of 1088 formal runo were conducted. Theae were fairly 
evenly divided between the various configurations and pilots, as indi¬ 
cated by Table 9. 

One of the first taskB of the simulation was to validate or modify 
the ratings scale and debriefing questionnaire to the satisfaction of 
the evaluation pilots. The final versions were used throughout the 
remainder of the simulation. Accordingly, this section starts with the 
rating acale assessment and then presents the airframe configuration 
assessments. 


127 






FOR PILOTED 





































TABLE 9. NUMBER OF RUNS 


CONFIGU¬ 

RATION 

PILOT 

TOTAL 

RC 

JF 

RH 

A 1 

59 

79 

79 

217 

a 2 

68 

48 

79 

195 

B 

61 

94 

58 

213 

C 1 

44 

77 

40 

161 

c 2 

53 

62 

58 

173 

D 

35 

62 

32 

129 

TOTAL 

320 

422 

346 

1088 


A. RATINGS SCALE AMD QUESTIONNAIRE ASSESSMENT 

The debriefing questionnaire was answered at the completion of the 
familiarization runs and again at the completion of the first set of 
tracking runs. During this period modifications were worked out for the 
rating scale* The modified scale was then used for the final Bet of 
tracking runs. 

1. Questionnaire 

The debriefing questionnaire proved essential to obtaining a consis¬ 
tent set of information for interpretation of results. However, the 
simple yes or no answers requested did not suffice. The tasks and 
resulting aircraft motions were so complex as to require considerable 
explanation. For example, more than a single mode of departure or spin 
was observed in some cases. Departure warning and severity also varied 
with rate of stall onset (&) and pilot usage of controls. Thus accom¬ 
plishment of the questionnaire often became quite time-coneuming, with 
several pages of accompanying explanation. This was augmented by tape 
recordings made, on-line, of pilot commentary while performing the simu¬ 
lation runs. As a consequence, pilot responses were so lengthy — and 







in so many cases repetitive - that it is not practical to incorporate 

them in this report. 

However, the basic conclusion reached by both pilots and experi- 
mentors was that such a questionnaire is highly desirable for stall/ 
departure/recovery investigation. The only change in format deemed 
necessary is provision of more space to accommodate lengthy answers. 

2. Rating Scales 

As in the previous simulation (Reference 5), the pilots felt that 
the Cooper-Harper rating scale was inappropriate for stall/departure/ 
recovery evaluation where aircraft characteristics rather than task 
performance are being assessed. On the other hand, the overall concept 
of the Figure 72 departure ratings scale was appreciated, although somu 
alterations were required before this scale met with full approval. In 
particular, the pilots did not like the open-ended aspects of the scales 
(the scale extends beyond the first and last descriptors). This clearly 
gave problems, and they preferred that the scale be bounded with a 
simple descriptor at each end. Second, the time required to recover and 
the control application necessary to effect recovery needed to be sepa¬ 
rated. Recovery controls also deserved a further breakdown, with one 
rating for required control complexity and a second rating for control 
application timing. Finally, it was decided that an assessment of 
hazard was difficult to make because this changes with situations such 
as low versus high altitude. The pilots preferred to provide an overall 
assessment baaed upon the acceptability of the flying characteristics at 
high angles of attack. The end point descriptors selected by the pilots 
were, in their jargon, "Sierra Hotel," meaning they would like to fly 
the aircraft at high AOA and "Delta Sierra," meaning unacceptable flying 
characteristics. The resulting rating scale is shown in Figure 73. 
Both test pilots were happy with this final form and considered that it 
covered all key factors, minimized the descriptor conflicts, and did not 
constrain the pilot's quantitative evaluation. 

The rating scale of Figure 73 was employed in the final set of 
tracking tasks. Resulting characteristic assessments for the six 
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vehicle configurations are summarized In Figure 74. The numerical 
racings shown for each pilot are the averages of separate ratings given 
for each of the three tracking tasks (Tj_ through T3). Considering the 
small statistical run sampling and the fact that all configurations were 
intended to have basic handling deficiencies, the results are considered 
to be very gocd. It is to be expected that all ratings would tend 
toward the poor (high numerical value) end of the scales and they gener¬ 
ally are. Most importantly, they also reflect the specific dynamic and 
handling variations intended to be exhibited between configurations and, 
in most instances, the ratings of RC and JF demonstrate agreement in 
variational differences if not in absolute levels. The ratings of pilot 
RH tend to vacillate between those ot RC and JF and also (except for 
recovery control timing) tend to agii.ua with the poorer rating. This 
downrating probably reflects pilot RH's relative inexperience in depar¬ 
ture/ ipin/recovery and the violent PSGs which high-performance fighters 
often exhibit. Thus the ratings of RC and JF are considered the more 
significant. 

Tho separation between RC and JF ratings bn warning and motion 
severity is consistent with differences in piloting techniques. Gener¬ 
ally JF is about one rating point harsher than KC because of a more 
rapid pull into departure and spin. This will be discussed in more 
detail later. 

Recovery control timing assessment shows very good agreement between 
the pilots except for Configuration C2. This primarily reflects timing 
of control release . The differential is about one rating point or less 
except Configuration C£ where RC had difficulty with the two different 
spin modes the aircraft could have. Timing of control release deter¬ 
mined whether an oscillatory (recoverable) or flat (unrecoverable) spin 
was obtained. 

Time to recover and overall rating show almost Identical trends and 
thus reflect predominant concern for loss of altitude in any departure/ 
spin situation. Configuration C^ had an oscillatory, recoverable 
spin. Configuration D did not spin but could end up with appreciable 
altitude loss if the pilot kept fighting the "wallowing" departure mode 
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(described in Subsection II.B). All other configurations had flat spin 
modes that were difficult if not impossible to stop. The significant 
disapproval of the Configuration D pitch-up is also evident for pilot RC 
on both of these ratings. 

The rating scale for recovery control complexity may need additional 
refinement. Pilot JF rated all configurations at 1 because he merely 
released controls to initiate recovery and. by definition, this is as 
simple as things can be. Pilot RC, with less severe departures, was 
able to discern some influence of the Btick and rudder in effecting 
recovery and thio was also reflected in his being more critical of 
recovery control timing. 

However, aspects pertaining to recovery (control usage and recovery 
time) for configurations ocher than D are not necessarily significant 
because the simulation cannot be considered valid for spin and spin 
recovery. The motions were generally so violent and the out-the- 

windscreen sky/earth display so featureless that the pilots had to 
resort to the all-attitude ball for recovery. This is not completely 
unnatural for the military pilots, but did make the task more diffi¬ 
cult. Additionally, the aerodynamic coefficients employed in the 
simulation do not represent steady spin conditions. Probably the most 
that can be attributed to recovery evaluations is that they exercised 
the rating scales. 

In summary, the departure rating scales of Figure 73 were accepted 
and supported by the pilotB and, based upon this small statistical 

sample, the numerical ratings given reflect characteristics "designed 
into" the vehicle configurations. The spread in ratings between the Air 
Force test pilots (generally about O.S on a 5-point scale) is consistent 
with spreads normally expected in Cooper-Harper ratings (about 1.0 on a 
10-polnt scale). Greater spreads generally reflect sensitivity to 
piloting technique; that sensitivity itself may be an attribute of the 
scales. 

Since vehicle attributes rather than a task performance level 

are being rated, use of the rating scales must be accompanied by a 
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qualitative assessment/description. The qualitative information should 
include, as a minimum, the following: 


• Warning 

- Type 

- Clarity 

- Margin 

• Departure 

« Resistance (susceptibility) 

- Type 

- Severity 

Ability of pilot to delay or prevent 
• Control action taken 

e Demands on the pilot 

• Post-Departure Motion 

- Type of aircraft motion 

- Severity 

• Recovery 

- Rapidity 

- Recovery controls 

- Demands 

« Ability to recognise 

a Ability to perform necessary control action 

B. COMFIGUBATIOM ASSESSMENTS 

As noted previously, the goals were to validate the influence of the 
key aerodynamic coefficients in determining departure characteristics; 
evaluate the influence of varied maneuver-limiting factors on high AOA 
maneuvering control; and identify potential flying quality criteria 
in terms of departure resistance, warning, and severity for possible 
application to the flying quality specification. In this subsection, 
predicted vs. observed characteristics aie reviewed and demonstrate 
that the above goals have been achieved, although therw were seme sur¬ 
prises in the process. The unaugmented airframe configurations will be 
assessed firot, and then the augmented cases. 
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It is necessary to rely heavily on qualitative statements contained 
in the debriefing questionnaire and commentary taped on-line as the 
pilots talked themselves through the runs. The questionnaires were 
filled in only during the initial progression through the test matrix. 
At that time the pilots were "feeling out” the various configurations 
and produced voluminous commentary on each, covering the primary obser¬ 
vations plus variations due to different control applications, tech¬ 
niques, rates of onset, tracking maneuvers, etc. The assessments 
reported in the following paragraphs have been distilled to the simplest 
possible factors consistent with the goals of the simulation. 

Only assessments of the fighter test pilots are reported because RC 
and JF are most experienced in stall/departure/spln testing and were 
more observant of differences between configurations. It became 
apparent early in the data analysis that the two test pilots were 
employing almost opposite approaches to stall/departure. RC employed a 
cautious, slow Increase in AOA and generally was able to detect subtle 
changes in vehicle stability or response characteristics. JF used an 
aggressive, rapid AOA increase which did not give time to detect such 
warnings. As a consequence, he generally pulled to a higher AOA with 
leas speed bleedoff and obtained more severe post-stall gyrations 
(PSG). Thus, differences in pilot closed-loop control techniques pro¬ 
duced a confounding influence which must be recognized before delving 
into other details! departure warning and severity are a function of 
pull-up rate and lateral-directional controls application immediately 
preceding departure. 

1. Unang—ated Airframe 

a. Configuration Aj 


Predicted versus observed pre- and post-departure characteristics 
are shown in Table 10. The warning Pilot RC observed with increasing 
AOA Is consistent with the graceful (gradual) degradation suggested 
by the open-loop dutch roll and roll numerator root migrations of 
Figure 55. It appears the negative C m ^ contribution was observed as a 
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TABLE 10. CONFIGURATION A^ MANEUVER-LIMITING FACTORS 


PREDICTED 

OBSERVED 

RC 

JF 

Warning: 

Wing rock 

(16 < o < 22 deg) 

Wing rock low ampli¬ 
tude, high frequency 

None 

Roll reversal 
(ct > IB dag) 

Roll revcraal 

(a > 20 deg); g-break 

(20-22 deg a) 

None; g-break 


Wing rock large ampli¬ 
tude, low frequency 


Rating: 

2.8 

4.8 

Departure: 

Noaa allca 
(a > 22 deg) 

. 

Noae allca followed by 
roll 

Nose slice and roll 

Poet etall gyratlona 

Large yaw, pitch, roll 
oaclllation and apln 

Violent ooclllatlona 
and spin 

Severity Rating: 

2.2 

4.6 

Departure/Spin 
Susceptibility 
(Welsaman criterion): 

o < 21 deg: mild or 
no departure 

Release a * 20 deg; 
yaw initially alow, 
builds rapidly 


21 < a < 24 deg: 
moderate rolling 
departure 


Release a • 23 deg; 
abrupt nose slice, 
bullda exponentially 

a > 24 deg: atrong 
rolling departure, 
high apln susceptl- 
bility 



Rating: 

Reaiatant 

Extremely Susceptible 

Commenta: 

Maneuvers tend to mask 
warning; not masked if 
pay attention 

Most violent poat 
stall gyration oscil¬ 
lations of any con¬ 
figuration 
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"g-break." The nose slice departure is as predicted, but the onset is 

noticed prior to C nft becoming negative. Thus onset might be influ- 
p dyn 

enced by pilot control, and thus associated with negative l/T^. Pilot 
RC considered the configuration departure resistant (R) because of ample 
warning and his tendency to initiate recovery at a * 20 deg where the 
Welseman criterion predicts mild or no departure. Pilot RC therefore 
gave fairly good high AOA ratings. 

Pilot JF eaw no warning except the "g-break." He obtained an abrupt 
noBe slice and violent post-stall gyration (PSG). The tendency to flat 
spin may be due to the large negative C n ^ and small Cat hia recovery 
initiation point and the further destabilizing influence of sideslip on 
the denominator and numerator roots, all of which portend high body-axis 
yaw rate and little roll. Inertia effects would tend to raise the nose 

and align the yuw rotational axis with the velocity vector-a flat 

spin. 

This configuration also exhibited violent oscillatory spin charac¬ 
teristics which could be reduced if a large nose-down pitch rate devel¬ 
oped as the nose was slicing. This could be produced via forward stick 
coinciding with the "g-break." The resulting rq Inertia cross-coupling 
produced a large roll acceleration and quite wild gyrations (s.g., 
Figure 75). Pitot JF noted this configuration to have the most violent 
PSG characteristics of any configuration and rated it extremely suscep¬ 
tible (ES) to departure. 

While the Weissman criterion correctly predicted the spin suscep¬ 
tibility, the initial departure was in yaw rather than roll. The 
maneuver-limiting factors are the strong nose slice and Bpin character¬ 
istics. The causal factors are combined open- and closed-loop (roll) 
instabilities accompanied by signifeant negative C^. 

b. Configuration B 

As predicted, this configuration is similar to Configuration A, but 
with less departure warning (aee Table li). The low-frequency yaw/roll 
oscillation that Pilot RC observed above 21 deg AOA is consistent with 
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TABLE 11. CONFIGURATION B MANEUVER-LIMITING FACTORS 


PREDICTED 

OBSERVED 

RC 

JF 

Warning: 

Roll ravaraal 
(o > 18 dag) 

Resists roll command; 
g-break 

Low frequency yaw/roll 
21 deg a 

No warning 

Rating: 

4.2 

4.8 

Departure: 

Noaa allca 
(a > 22 dag) 

Poat atall gyratlona 

Nose slice 

Violent roll reversal 
or spin 

Nose slice 

Severe and violent 
oscillation 

Severity Rating: 

2.7 

4.5 

Departure/Spin 
Susceptibility 
(Weisaman Criterion) 

a < 21 deg: mild or 
no departure 

Release a * 21 dag; 
yaw followed by roll 
and pitch down 


21 < a < 24 deg: 
moderate rolling 
departure 

a > 24 deg: strong 
noae slice, high apln 
susceptibility 

-- 

Release o * 24 deg 

Nose drift followed by 
roll; prone to flat 
spin 

Rating: 

Susceptible 

Extremely Susceptible 

Comments: 

Cannot discern onset 
of departure; Insidi¬ 
ous 

Low acceleration in 
noao slice 
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open-loop dutch roll dynamics (see Figure 57). The lower frequency 
resulted In less acceleration Into the nose slice. Coupled with roll 
reversal, this apparently masks the departure onset. Although he Ini¬ 
tiated recovery at about the same AOA and thought the departure severity 
(pilot rating ■ 2.7) Is about the same as Configuration Aj, he consid¬ 
ered Configuration B to be departure susceptible (S) because of the poor 
warning (PR ■ 4.2). Pilot JF could see little difference between A^ and 
B except the increased tendency to flat spin due to the combined high 
aerodynamic roll rate damping and strong directional divergence. He 
again rated it extremely susceptible (ES). 

This configuration exhibited all anticipated flying characteris¬ 
tics. The Welssman criterion correctly Identified departure and spin 
severity but not type. Again, the maneuver-limiting factors are the 
severe nose slice and flat spin tendencies* The causal factors are the 
open- and closed-loop directional Instabilities aggravated by high 
negative C^. 

c. Configuration 

Configuration (Table 12) with positive dynamic stability (due to 
large negative Cj^) had a much greater tendency to develop roll oscilla¬ 
tions. Both pilots noted that any lateral control input excited wing 
rock and produced a tendency to lateral PIO. Since this sensitivity 
persisted over a wide AOA range, it was not considered a warning but 
actually a masking of departure onset. Both pilots considered the 
warning to be poor (PR > 4). 

This configuration is stable open-loop and driven unstable and into 
departure by closing the attitude loops. Departure characteristics are 
less consistent (predictable) because the nature of departure depends 
upon pilot control activity* It exhibited at least two and possibly 
thrse different types of departure onset, which further aggravated the 
poor warning. The lateral phugoid mode could contribute the additional 
departure mode (see Figure 59). One pilot detected more tendency to yaw 
and the other more tendency to roll. However, they generally were in 
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TABLE 12. CONFIGURATION C 1 MANEUVER-LIMITING FACTORS 


PREDICTED 

OBSERVED 

RC 

JF 

Warning: 

Wing rock 

(16 < a < 26 dag) 

Wing rock: any 
disturbance 

Wing rock: tendency 
to low frequency 
pilot-induced 
oscillation 

Roll reversal 
(a > 18 deg) 

Low frequency, high 
amplitude oscillation 
masks departure 


Rating: 

4.3 

4.0 

Departure: 

Roll departure 

Post-stall gyrations 

Nose slice 

Roll 

Two distinct types: 

a) roll, little yawi 

b) yaw, followed by 
roll 

Severity Rating: 

1.6 

3.7 

Daparture/Spin 

Susceptibility: 

(Weissman Criterion) 

a < 21 deg: no depar¬ 
ture 



a > 21 deg: mild 
rolling departure, low 
spin susceptibility 

Release a * 32.5 deg, 
nose slice and roll 

Release a * 25.5 deg, 
milder departure, 
oscillatory spin 

Rating: 

Susceptible 

Susceptible 

Comments: 

Takes tine to see yaw 
ie not going to stop 

Not prone to flat spin 
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closer agreement for this configuration concerning warning, nature, and 
severity of and susceptibility to departure, probably due to the posi¬ 
tive dynamic stsbility which prolonged the time period the airplane 
could be maintained under semi-control. Thus, both pilots had suffi¬ 
cient time to observe its high AOA characteristics fully before closed- 
loop instability predominated. 

Both pilots considered the departure and FSG severity to be somewhat 
lees than for either A^ or B. The configuration was spin-prone, but 
spins were generally highly ocillatory with some possibility of recov¬ 
ery. It was not prone to flat spins. Departure and spin characteris¬ 
tics and susceptibility are in agreement with the Weissman criterion. 

The initial or basic maneuver-limiting factor is the strong wing 
rock tendency which extends from relatively low AOA through departure 
and is aggravated by closed-loop instability (roll PIQ). Departure is a 
secondary factor. Warning may also be adversely influenced by the low- 
frequency, large-amplitude, lateral phugoid oscillation. 


d. Configuration D 


This configuration (Table 13) had less directional instability but 
did have sideslip-induced pitch-up. The expected wing-rock tendency 
from a lightly damped dutch roll was not observed by either pilot. Both 
noted departure onset as a slow, persistent yaw (noBe slice) accompanied 
by mild roll reversal end pitch-up. The onset wee sufficiently mild 
that departure wae eemi-contTollable throughout the achievable AOA 
range. A. poor departure warning rating (PR > 3.5) resulted. 


Open- and cloead-loop instability parameters (Wj and u)2 or C r 


▼ H ayn 

and LCDP) are small for this configuration and lead to sloppy control 
and low-frequency wallowing which masks daparturs. At the same time, 
tha wallowing does not generate sufficiently rapid motion to excite 
inertia croee-coupllng and PSG. All pilots tended to continue fighting 
to maintain control well past full stall, incurring excessive altitude 
loss. However, if controls wars released at any time the aircraft would 
immediately go into a noae-low spiral and recover by itself. 
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TABLE 13. CONFIGURATION D MANEUVER-LIMITING FACTORS 


PREDICTED 

OBSERVED 

RC 

JF 

Warning: 

Wing rock 

<16 < a < 26 deg) 

Roll reversal 
(a > 20 deg) 

Yaw, roll reversal, 
pitch up 

Sloppy control, large 

6 and $ overahoots 
mask departure 

Large definite yaw 
followed by reversal 
and pitch up 

Rating: 

4.0 

3.5 

Departure: 

Pitch up 

Poet stall gyration 

Nose slice/pitch up 

None, never spin 

Nose slice/pitch up 

Mild, no spin 

Departure/Spin 
Susceptibility 
(Welssman Criterion): 

All a: no departure 

Release a * 27.6 deg, 
slow nose slice, pitch 
up, some roll 

Release a * 27.6 deg; 
slow consistent nose 
slice, roll, departure 
nose high 

Rating: 

Resistant 

Resistant 

Comments: 

Do not like pitch up 
at departure; cannot 
get nose down 

Always recoverable if 
release controls 






























Departure warning was considered poor (e.g., 4.0) by RC because of 
the wallowing, but at the sane tine the configuration was considered 
departure resistant (R)» Pilot JF assessed It as resistant primarily 
because it would not spin. Pilot RC's principal objection to this 
configuration was the positive C^. Even though it eliminated the 
severe pitch-down in nose slice and the resulting inertia-coupled PSG 
exhibited by other configurations, ha did not like the nose rising, no 
matter how mildly, when in or near departure conditions. 

The maneuver-limiting factor for this configuration is the wallow¬ 
ing - a general lnsansitivity to or sloppiness of response to pilot 

control inputs prior to the departure AOA. The key causal factors are 
considered to be the relatively small (due to small % and C&g) 
which produces very low DC gain in the roll frequency (Bode) domain and 
positive which apparently contributed to overcontrol in pitch. 

e. Summary 

The intended high AOA departure characteristics were obtained in the 

four unaugmented airframe configurations. The results demonstrated the 

Influences that changes in the six key aerodynamic coefficients c&n have 

on departure warning, severity, and susceptibility. The .su’csptibility 

of Configuration C^ to departure was somewhat revealing: the open-loop 

parameter C n . , by Itself, is not a sufficient indicator of departure 
s dyn 

resistance. However, when combined with the closed-loop parameter LCDP 
(or «^), the two are quite accurate indicators of departure suscepti¬ 
bility and severity as advanced by Welssman. 

Departure severity also appears influenced by the static cross¬ 
coupling derivatives <£&> Na» and Mg. Severe departure and PSG were 
obtained for the two configurations in which combinations of these 
derivatives shifted lateral-directional denominator and numerator roots 
into the right half-plane of the root locus plot with S i< 0 (Figures 55 
ai;d 57). Leso severe departure was obtained for combinations which 
primarily shifted numerator roots info the right half-plane (Figure 59), 
and the mildest departure was obtained for combinations which produced 
very little shift of lateral denominator and numerator roots with 
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sideslip (Figure 61). The literal approximate factors developed in 
Section III indicate that the small shifting is due to the combined 
Influence of small N£ and positive However, the nose rise asso¬ 

ciated with large positive Mg was disconcerting to one pilot. 

The simulation results also demonstrated that wing rock or lateral 
PIO is not a satisfactory warning of impending departure if ''.here is 
sufficient elevator power to pull through the warning region rapidly or 
if the AOA range in which it occurs is too wide. 

2. Augmented Airframe 

a. Configuration A 2 


The flight control augmentation was configured to provide a high 
level of roll rate damping, and to improve roll maneuvering response via 
roll rate command (p c ) and reduced adverse aileron yaw (SRI). This 
augmentation was also expected to eliminate all cues of impending depar¬ 
ture. The departure was anticipated to be a strong nose slice due to 
increasingly negative C n g as AOA increased and, once sideslip started to 
build, it was expected that the characteristics would be very similar to 
those of Configuration 

Results summarized in Table 14 for Pilot JF are very much as antici¬ 
pated. However, Pilot RC detected a pre-nose-sllce resistance to roll 
commands (at AOA above approximately 21 deg) when attempting to track 
the maneuvering target. This resistance was observed only during the 
initial tracking sets and resulted in the pilot's neutralizing all con¬ 
trols before the directional Instability became severe (rating 1.8). 
The aecond set of tracking runs were done more hurriedly, the warning 
was not obtained (rating 4.3), and more severe departures ensued (rating 
3.7). The susceptibility assessment as resistant (R) by Pilot RC was 
given only on the initial set of runs. 


vi 


w 


The maneuver-limiting factor for this configuration Is the nose 
slice. But. since there la no warning, an AOA limit wculd have to be 
eet to provide a safety margin against inadvertent departure. The key 






TABLE 14. CONFIGURATION A 2 MANEUVER-LIMITING FACTORS 


PREDICTED 

OBSERVED 

RC 

JF 

Warning: 

None 

Resists above 

21 deg a 

Cluar warning with 
target 

None 

Rating: 

2.8 

(4.3 on reruns) 

4.8 

Departure: 

Nose slice 

Nose slice 

Nose Blice 

Roll and severe PSG 

Severity Rating 

1.8 

(3.7 on reruns) 

4.5 

Departure/Spin 
Susceptibility 

a < 23 deg: no depar¬ 
ture 



a > 23 deg: nose 
slice and revert to 
Configuration A| 
characteristics 

Release a * 25 deg; 
smooth, slow yaw 
followed by roll 

Release a A 32 deg; 
exponential yaw diver¬ 
gence tendency to flat 
spin 

Retina: 

Resistant 

Extremely Susceptible 

Comments: 

Depart without pilot 
aware; only cue is 
resistance to $ c 
against target 

Easier to fly 

Rate of yaw dependent 
upon 4 or g 
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causal factors are the strong static directional divergence and the 
elimination of departure warning by the augmentation system. 

Time traces of a Configuration A 2 departure from a winga-level pull- 
up are shown In Figure 76. The lateral Btick trace reflects about the 
same activity prior to departure as in Figure 75 but far leas sideslip 
and roll oscillation. After pilot neutralization of controls, the aug¬ 
mentation loops continue to deflect the lateral and directional control 
surfaces and apparently aid recovery (again compare with Figure 75). 
Interestingly, the actual departure and PSG motions of Figure 76 are 
remarkably similar to Reference 38 flight traces for the F-4E with 
augmentation systems on (Figure 77). Prior to departure the traces of 
the two figures are not comparable because the stick to rudder intercon¬ 
nect in the simulation reduced sldoalip excitation and provided better 
roll attitude control. 


b. Configuration C 2 


Based on the positive C n ^ and Improved (or LCDP), this aug¬ 
mented configuration was expected to have the best high AOA flight 
characteristics and little or no departure tendency. Xt did allow 
both pilots to consistently maintain control to higher AOA than did any 
other configuration, though abrupt nose slice departures were common 
(Table 15). In this case the roll damper countered the dynamic 
stability contribution and allowed the large negative C n ^ otatlc insta¬ 
bility to predominate. The aircraft responses by themselves provided no 
warning of Impending departure, although Pilot RC was able to extract 
some low-frequency PIO cues with respect to the target aircraft. He 
rated the configuration departure resistant (R) because of this cue and 
the mild (to him) nose olice motion. Without ths target he rated the 
configuration departure eusceptible (S), as did Pilot JF. 


This configuration was also spin-prone and, similar to C^, exhibited 
two dletlnct spin modes as noted in Table 15. It appeared to be more 
prone to spin flat than was the unaugmented configuration (C^), because 
the roll rate CAS automatically opposed any uncoramandsd roll. 


148 






















































TABLE 15. CONFIGURATION C 2 MANEUVER-LIMITING FACTORS 


PliEDICTED 

OBSERVED 

RC 

JF 

Warning: 

None 

Low frequency pilot- 
induced oscillation; 
g-break 

None, very steady 

Rating: 

3.2 

4.8 

Departure: 

None 

Nose slice followed by 
roll and pitch down 

Abrupt nose slice 

Severity Rating: 

1.5 

3.8 

Departure/Spin 

Susceptibility: 

Reslitant 

Release a * 31 deg; 
slow nose slice; 
almost always spin 

Release o * 30 deg; 
abrupt nose slice; two 
Bpin modes: 

a) oscillatory 
(recoverable) 

b) flat (nonrecover- 
able) 

Rating: 

Resistant 

Sueceptible 

Comments: 

Resistant because tar¬ 
get gives cue; without 
target is Susceptible 

Better flying quali¬ 
ties prior to depar¬ 
ture 











The maneuver-limiting factor for this configuration is the nose 

slice departure. The positive C na and ui^ throughout the achievable 

Pdyn t 

AOA range did allow it to be controlled to considerably higher ADA than 
the other configurations. In this respect it was less susceptible to 
departure. However, at the higher AOA directional instability even¬ 
tually required full opposing rudder deflection, which was immediately 
followed by departure. 


c. Overall Assessment 

The augmented flight control system designed to enhance character¬ 
istics in normal flight improved pre-departure control precision, 

reduced departure warning, and increased the tendency to flat spin. The 
roll rate CAS and the SRI together improved lateral-directional control, 
which allowed higher AOA to be reached before sideslip became signifi¬ 
cant. The major portion of this benefit derived from the SRI. For 

Configuration C 2 the roll rate damper suppressed the <Cg5 contribution 
to dynamic stability and thus increased the nose-slice departure and 
flat spin tendencies. Thus, this augmentation function actually 
degraded high AOA characteristics. 

C. SUMMARY 

The intended high AOA departure characteristics were obtained in the 
four bare-alrframe and two augmented configurations. The results demon¬ 
strate the strong Influence changes in the six key aerodynamic coeffi¬ 
cients (£q» f 3> «£p> Na> «nd Mg) can have on departure warning, 

severity, and susceptibility. The simulation plus literal approximate 
factors of Section III show: 

• The most violent departure characteristics and 
flat spin tender Lea to be associated with large 
negative N^, Ng, and Mg. 

• The least violent departure characteristics and 
greatest spin resistance to be associated with 
•mall to moderate negative and Kg and positive 
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• The tendency tc Inconsistent departure/spin char¬ 
acteristics (possibly several modes) to he asso¬ 
ciated with positive C n @cly n brought about by 
large ..negative Co Q and continued large adverse 
y.w <«$. LCDP )4 H 

However, the observed high AOA stall, departure, and spin dynamic char¬ 
acteristics also depended significantly on pilot control technique. 
This is consistent with actual flight experience in which it haB been 
observed, for example, that spin characteristics and recovery techniques 
can be quite different, depending upon the maneuver from which the spin 
occurred. 

The final high AOA departure/spin rating scale of Figure 73 was 
accepted and supported by the pilots. Results of this first application 
showed consistent reflection of both the airframe dynamic characteris¬ 
tics and any sensitivity to difference in piloting technique. However, 
the ratings must be accompanied by detailed qualitative description of 
the characteristics observed, i.e., what is being rated. 

In this vein, it might be noted that there are some differences 
between ratings for the Initial familiarisation and tracking runs 
reported in Subsection B and for the final tracking runs reported in 
Subsection A. These are due to several factors. First, Pilot RC was 
working with the original open-ended scale (Figure 72) during the ini¬ 
tial runs and the cloaed-ended scale (Figure 73) during the reruns; 
Pilot JF used the close-ended scale for both sets of runs. Second, 
during the first sequence of runs Pilot RC was unhurried in investi¬ 
gating the characteristics of each configuration. However, for the 
repeat runs he was under some pressure to complete the series of maneu¬ 
vers and configurations within time constraints. Thus, the rate at 
which he approached stall/departure was more rapid and this influenced 
his assessment of warning and motion severity, as demonstrated in 
Table 16. The greatest differences are for the augmented flight control 
configurations, where RC saw considerably less warning and more severe 
departure during the reruns. Results for JF, shown for comparison, 
reflect quite consistent ratings between the sets, since he always used 
a rapid pullup technique. 
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TABLE 16 


COMPARISON OF PILOT RATINGS, 
INITIAL VS. RERUNS 


CONFIGU¬ 

RATION 

RC 

JF 

WARNING 

MOTION SEVERITY 

WARNING 

MOTION SEVERITY 

INITIAL 

RERUN 

INITIAL 

RERUN 

INITI4L 

RERUN 

INITIAL 

RERUN 

A 1 

2.8 

2.9 

2.2 

3.0 

4.8 

5.0 

4.6 

4.3 

B 

4.2 

3.9 

2.7 

3.1 

4.8 

4.6 

4.5 

4.8 

Cl 

4.3 

3.6 

1.6 

2.0 

4.0 

4.6 

3.7 

3.9 

D 

4.0 

3.8 

3.5 

3.0 

3.5 

3.7 

1.5 

2.9 

A 2 

2.8 

4.3 

1.8 

3.7 

4.8 

5.0 

4.6 

4.5 

C 2 

3.2 

4.2 

1.5 

2.7 

4.8 

4.1 

3.8 

3.9 
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sectiob vii 


DEPARTURE SUSGEPTIBLITT CKITK&IOB 


Pitch control techniques and aggressiveness had obvious influences 
on the departure warning and susceptibility assessments shown in the 
previous section. When asked about specific cues and personal criteria 
which would prompt the pilot to acknowledge departure and therefore 
to initiate recovery, vague answers were obtained, e.g., "persistent 
or continued uncoomanded roll or yaw" or "suddon, rapid uncommanded 
motion — somewhere around 20-30 deg/secHowever, detailed scrutiny 
of the time traces failed to reveal correlation between release of back 
stick and various motion quantities (p» r, 8, etc.). These quan¬ 

tities were always changing rapidly and timing became a big factor. 
Deviations of fractions of a second in initiating recovery (or reading 
the tracao) produced a large variation in the motion amplitudes (e.g., 
see Figures 57 and 58). The greatest correlation was found in AOA 
itself — posaibly because its rate of change was the slowest of all and 
therefore could be read and correlated with the greatest accuracy with 
what the pilots saw. 

Since the piloting task required closed-loop attitude control,, and 
there was considerable evidence of LCDP or wj being a key parameter,, 
attention was turned to the possibility of a closed-loop criterion. A 
simple analysis of theoretical roll loop stability was found to agree 
quite well with frequencies, damping ratios, and divergences obtained 
from the time traces. Furthermore, the results corroborated the differ¬ 
ences in pilot techniques and obaerved warning reported in the previous 
section. 

This section summarises the relationship obtained between the 
closed-loop analysis and the piloted simulation time traces and presents 
a closed-loop criterion which is consistent with and explains the seem¬ 
ingly diverse departure susceptibility assessments given by the two Air 
Force test pilots. 
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A. CLOS&D-LOO? ANALYSIS 


The lateral stick time trace of Figure 57 (Configuration k±) Is 
typical of most runs for Pilot RC. The phasing between 6 excursions and 
stick deflections Indicates little, if any, lead being generated during 
the PIO at a - 16 deg and above* Therefore, the simple block diagram of 
Figure 78 was chosen to model the unaugmented airframe, control system, 
and pilot. A Bode-root locus closed-loop survey plot for Configuration 
Aj at o ■ 16 deg and 6-0 deg is presented in Figure 79. Closed-loop 
gains and frequencies obtained from strip chart recordings of simulation 
runs for Pilot RC are shown on the Bode amplitude plot by the several 
horlsontal lines (K^) and tic marks (-4-). [For example, at 16 deg AOA 
the traces of Figure 57 show a pilot gain * - 4.5 (13 dB) and fre¬ 
quency w * 1.8 rad/sec«] The frequency of the undamped oscillation 
obtained from the strip chart is in remarkable agreement with the fre¬ 
quency for 180 deg phase of the analytical predictions, validating the 
assumptions that this pilot is adopting neither leed nor lag in attempt¬ 
ing to control roll attitude. Rather, he is merely reacting with a 
normal time delay with stick proportional and opposite to bank angle 
excursions. 



Figure 78. Assumed Pilot Lateral Loop Closure; 4 ♦ 5*tk 
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The measured gains range from about 10 to 20 dB. This ^riatlon. 
Identified on the root locus plot of Figure 79, indicates but small 
Influence on dutch roll divergence. The 13 dB gain point indicates a 
mild divergence which is also in agreement with the traces of Figure 57 
at 16 deg AOA. 

As AOA is Increased in Figure 57, the lateral stick amplitude 
remains fairly constant but the 4> excursions increase. Thus, the loop 
gain (K^) decreases slightly, the oscillation period increases to about 
6 seconds (w * 1 rad/sec), and a yaw rate divergence develops. Again, a 
closed-loop survey plot for a ■ 23 deg (Figure 80) shows excellent 
agreement, i.e., uij * 1 rad/sec and a first'-order spiral divergence 
of T; * 2 sec which results from the spiral mode being driven toward the 
RHP xero l/T^. 

Similar trends were obvious from Pilot RC run traces for moat con¬ 
figurations. He appeared to be following the neutral stability boundary 
of the oscillatory mode and initiating recovery when the first-order 
divergence mode exceeded soma as yet unidentified divergence rate. 

Therefore, using the same pilot model, the maximum boundary for 

stable modes was calculated and plotted as shown in Figure 81. ThiB 

shown two boundaries. The upper is the gain not to be exceeded at 8 « 0 
and the lower is the reduced limit gain for 8 " 5.5 deg. To the right 
(higher AOA) of either boundary a lateral or directional instability 
exists and if AOA is increased K^, must be decreased to maintain sta¬ 

bility. 

The roll loop gains employed by both pilots were determined from the 
time traces and located on similar stability boundary plots drawn for 
each vehicle configuration. Figure 82 is the plot for Configuration A^ 
again, but now showing gain bands for. each pilot. The band for RC shows 
that he tended to keep hia gain within 10 dB above the maximum stable 
boundary and to decrease his gain as AOA was increased. On the other 
hand, JF started with a lower gain and tended to keep it constant during 
the pull-up. 
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Amplitude (dB) 


Effeet of Pilot Loterol Loop Closure 




Figure 80. Survey Plot for Soil Loop Closure; 
a 0 - 23 deg, B 0 - 0 deg 
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Figure 81. Maximum Pilot Gain for Stable $ ®stk Closure; 
Configuration Ai 



Figure 82. Roll Loop Closure Gains Employed with 
Configuration Aj 












Figure 33 Is the gain plot fot Configuration A 2 • Again JF used a 
lower, constant gain while RC tended to follow the stability boundary 
until the sudden drop. Figure 84 is the plot for Configuration C^, 
showing results similar to those for Configuration A^. This configura¬ 
tion is somewhat unusual In that a second stable region exists for 8 f 0 
at very high AOA. Figure 85 for Configuration C 2 is similar to Configu¬ 
ration A 2 except there is no dropoff in stability boundary at high 
AOA. The results for Configuration D (Figure 86) are consistent: Pilot 
RC followed the stability boundary and hence should detect departure 
onset; Pilot JF used a lower, constant gain which suddenly crosses the 
boundaries as AOA is increased. He therefore obtained little warning. 
The plot for Configuration B is presented in Figure 87. This vehicle is 
so well damped in roll that it was difficult to obtain the necessary 
data from time traces except for one run by RC. 

At the bottom of Figures 82-87 negative values of the open-loop roll 
numerator root 1/T^ are plotted versus a (at 8 ■ 0). Figure 85 does 
not show this extra plot because is positive to o > 35 deg for this 
configuration. Also identified on each plot is the average AOA at which 
recovery was Initiated by Pilots RC ( ^symbol) and JF (O). It will be 
noted that: 

• In moat instances a decrease in the boundary 
is related to l/T^ becoming negative (RHP zero). 

• In all but one case recovery is initiated by both 
pilots after exceeding the K, boundary and exper¬ 
iencing negative l/T^,, (the 9 only exception is 
Configuration Cn for which the numerator roots 
always lie in the LHP and there is no dropoff in 

e Pilot RC always employs higher roll loop gain, 
always follows the Ka boundary, and generally 
initiates recovery at 9 lower AOA. However, for 
Configurations and C 2 he did explore higher 
AOAs before recovery. 

• Pilot JF tended to adopt a lower roll loop gain 
and kept it constant during his more rapid 
(aggressive) pull-up. He generally saw little or 
no warning and penetrated further into instabili¬ 
ty regions before initiating recovery. 
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0 Pilot RC 
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• O Pilot JF 



Figure 83. Roll Loop Closure Gains Employed with 
Configuration A 2 






Figure 84, Roll Loop Closure Cains Employed with 
Configuration Cj 
















Figure 85. Roll Loop Closure GainB Employed with 
Configuration C 2 



Figure 86. Roll Loop Closure Cains Employed with 
Configuration D 
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Figure 87. Roll Loop Closure Gain Employed 
with Configuration B 


The foregoing observations appear to form the basis for each pilot's 
definition of departure susceptibility (cf. Section VI): 

RC: Departure resistant (R) if warning is clear and 
consistent 

Departure susceptible (S) if warning is unclear, 
masked, or inconsistent 


1 


% 

1 

i 


JF: Departure resistant (R) if always recovered 

Departure susceptible (S) if sometimes recovered 

Departure extremely susceptible (ES) If never .'j 

recovered 

Thus each pilot is viewing different aspects of the departure. Fortu- ^ 

i 

nately, this difference «erves to increase the information potential of • 

the simulation program. This is probably typical of flight test, and | 

indicative of why such stall/spin programs should be flown by more than '] 

one pilot. | 
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B. POSSIBLE DEPARTURE SOSCBPTI1ILITT CRITERIA 


A comparison between departure/spin susceptibility predicted by the 
Weissman criterion end the eseeseaents provided by the two pilots Is 
shown is Table 17. Assessments substantially in agreement with predic¬ 
tion are shown in boxes. Obviously the aggressive pilot observed the 
worst possible characteristics of each unaugmented configuration as 
predicted by the criterion. The less aggressive pilot experienced 
something quite different. As noted previously, the criterion does not 
lend Itself to prediction of the frequency-dependent augmented airframe 
chercterlstlcs, but these configurations were expected to be less sus¬ 
ceptible to departure. This influence was observed only by the lees 
aggressive pilot. 


TABLE 17 

COMPARISON OF PREDICTED AMD ACTUAL DEPARTURE 
SUSCEPTIBILITY ASSESSMENTS 


DEPARTURE SUSCEPTIBILITY 

CONFIGU¬ 

RATION 

PREDICTED 

RC 

JF 

A 1 

ES 

R 


ES 


B 

ES 

S 


ES 


Cl 

S 

0 

S-ES 

D 

R 

0 


0 


a 2 

NA 

R 


ES 


C 2 

NA 

R 


S 



R - Resistant ES * Extremely susceptible 

S ■ Susceptible NA ■ Not applicable 
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One Important difference between prediction and simulation was the 
nature of departure. Figure 63 predicts predominantly rolling departure 
with no indication of yaw departure. Our Configurations A, B, and D 
exhibited initial yaw excursions sometimes followed by roll. Configura¬ 
tion C exhibited two or even three different departure modes which were 
dependent upon control application at onset of departure. However, a 
rolling type motion did predominate. Thus, the Welssman criterion was 
not as useful in predicting the nature of departure observed by the 
pilot. 

The closed-loop analysis to Identify causal factors behind the 
widely differing pilot ratings for departure/spin susceptibility pro¬ 
duced a strong relationship between pilot ratings for both pilots and 
penetration of the roll numerator into the RHP. Figure 88 shows the 
value of the real part of the numerator root at the Instant the pilot 
decided he had departed, or was about to depart, and Initiated recovery. 
These values are plotted against the AOA at which recovery was started. 
(Note this is not the usual root locus ju> axis.) The points represent 
all six vehicle configurations as evaluated by both pilots. 

The division between departure-resistant (R) and susceptible (S) 
ratings is seen to lie at roughly -0.5 rad/sec. This corresponds to a 
time to double amplitude of approximately 1.4 sec. Zeros which lie to 
the left of this line apparently limit the first-order divergence to a 
rate slow enough for pilots to respond and recover. Zeros to the right 
of the line apparently allow divergence rates so fast that the pilotB 
cannot prevent departure. Again this is consistent with the previously 
noted definition of departure whereby both pilots indicated a threshold 
on rate of motion; however, the pilots were vague as to the value (e.g., 
"maybe 20 or 30 deg/sec"). One data point in Figure 88 violates the 
boundary. This is the augmented Configuration C 2 > In this case the SRI 
eliminated adverse yaw and thus made the vehicle more departure resis¬ 
tant (as viewed by RC). However, it could be departed; and then the 
augmentation produced pro-spin control. Pilot JF rated spin suscepti¬ 
bility (and recovery) and apparently rated this configuration accord¬ 
ingly. 
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Closed Loop Divergence Potential 

Figure 88. Departure Susceptibility Rating Versus 
Lateral Closed-Loop Divergence Potential 


The Implication is that if the combined aerodynamics and flight 

control system design Is such that 1/T^ never exceeds -0.5 throughout 

the achievable A.0A range, the airplane will be departure-resistant. It 

should be noted that this criterion places no restriction on open-loop 

stability. For example, C n - can be negative and, In fact, is nega- 

Pdyn 

tlve for Configuration D (see Figure 62) which is rated departure- 
resistant (see Table 13) by both pilots. 

o 

Since < 0 ^ is the dimensional form of! LCDP end since, in general, 
ll/T^jl 2 ImJI, than one can relate the above 1/T^ boundary to an 
equivalent LCDP. For the flight conditions, Inertias, etc., employed 
in this simulation, 1/T<^ of -0.5 corresponds to LCDP of -0.001. This 
coincides with Welseaan'e boundary between Regions K end B for positive 
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C no , see Figure 89; however, it is a little more conservative at 
p dyn 

negative C n ^ . Thus, the results of our simulation are compatible 
with and support Weissman"s empirically derived LCDP boundary. The key 
difference in the criterion is that 1/T^ is not restricted to airframe 
lateral-directional static coefficients but can be applied to the com¬ 
pletely coupled 6 DOF airframe with a full complement of augmentation, 
stlck-to-rudder crossfeed, etc., throughout the aircraft development 
cycle. 


In addition, results of our simulation tend to indicate that another 
boundary might be appropriate in the upper left quadrant of Figure 89 to 
create two additional regions, E and F, as suggested in Reference 40. 
Region E would be classified as mild directional divergence and moderate 
spin tendency. Region F might be classified bb severe directional 
divergence and strong spin tendency. Note from Figure 63 that our air¬ 
craft Configuration D, which had mild departure characteristics but 
little or no spin tendency, would extend into Region E. 


LCDP 



Figure 89. Possible Modifications to Welsaman Criterion 
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Additional support for the limit on 1/T^ can be drawn from the 
current MIL-F-8785C spiral divergence limit* The Level 3,* Category A 
and C Flight Phases requirement is that time to double amplitude be not 
lees than 4 sec* This limit came, in part, from the Reference 41 inves- 
tigation in which aircraft configurations having various spiral diver¬ 
gence races were rated by twelve pilots in a landing approach task* 
Results of that investigation, shown in Figure 90, indicate that diver¬ 
gences having a time to double amplitude of less than 2.4 sec were 
considered intolerable. The time to double amplitude (T 2 ) can be 
converted to an equivalent first-order time constant (T) by the rela¬ 
tionship 


T 2 * 0.74T 

A plot of pilot rating versus spiral-mode inverse time constant (1/T) is 
presented in Figure 91. The current -3785C Level 3 spiral limit is 
shown along with the proposed 1/T^ limit. 

The two limits are in good agreement when it is remembered that 
1/T a represents a pole (modal response) while 1/T^j represents a zero. 
Under closed-loop roll control the migration of an airframe pole to¬ 
ward l/T^ would be dependent upon pilot gain; infinite loop gain would 
be required to obtain a closed-loop pole exactly at the zero location. 
A more realistic situation is shown in Figure 92, which is a Bode-aiggle 
plot for Configuration A^ at a ■ 21 deg, 8*0 deg. The dashed line 
terminating at l/T^ ■ -0.6 is the locus of the closed-loop first-order 
divergence root for increasing gain, . Two typical gain lines for 
Pilot RC extracted from simulation traces at a - 21 ± 1 deg are shown. 
These closures result in closed-loop roots at -0.2 and -0.33 rad/aec. 
Referring back to Figure 91, it may be observed that these straddle the 


A 

Level 3: “Flying qualities such that the airplane can be con¬ 

trolled safely, but pilot workload is excessive or mission effectiveness 
is inadequate, or both. Category A Flight Phases can be terminated 
safely, and Category B and C Flight Phcaas can be completed.” 
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Time to Double Amplitude of Aperodlc Motion, T.a (sec) 

Figure 90. Variation of Pilot RatingB with Time to Double Amplitude 
of Aperiodic Motion; Landing Approach Condition (from Ref. 40) 


-8785C 



Tolerable 


Intolerable 


Figure 91. Variation of Pilot RatingB with Equivalent 
Aperiodic Inverae Time Constants 
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1/T g boundary and are consistent with the Reference 40 division between 
tolerable and intolerable handlings Again, our Pilot RC considered 
Configuration A^ to be departure-resistant (R). 

The above discussion tends to show that the 1/T^ limit selected 
(-0.5) may not be conservative, since a very tight roll loop closure 
could result In closed-loop roots in the intolerable range of Figure 
91. Further investigation will be required to determine if the limit 
should be reduced. 

C. SUMMARY 


Analysis and piloted simulation have Bhown that: 

• The pilot's perception of departure suscepti¬ 
bility was found to be correlated with movement 
of one root of the roll numerator for lateral 
stick control into the RHP of the root locus, 
i.e., a nonminiraum-phase zero. If pilots could 
fly the aircraft to such AOAs, they rated high 
AOA configurations which produced a zero, 1/Tx , 
more negative than -0.5 departure-susceptible. 

If this boundary was not or could not be ex¬ 
ceeded, the aircraft was considered departure** 
resistant. This rating is not a function of the 
sign or magnitude of the dynamic stability param- 

Cter Sdyn* 

• A value of 1/T*. ■ -0.5 corresponds for the air¬ 
frame tested ter an effective LCDP of -0.001 and 
thus is consistent with and supports the empiri¬ 
cally derived LCDP departure boundary developed 
by Wetasman. However, results of the simulation 
were not in agreement with the types of departure 
predicted by Weiseman in that negative C n g d re " 
gions produced yaw or nose-Blice type departures. 

• A value of 1 /Ta. - -0.5 also is compatible with 
the current MIU-T i '-8785C Level 3 limit for spiral 
divergence iu Category A and C Flight Phases. 

It is suggested that as the aircraft design/development cycle pro¬ 
gresses past the static wind tunnel phasr the Welssman criterion for 
departure/spin susceptibility be replaced by a negative 1/T^ limit. 


172 






SECTION VIII 


CONCLUSIONS AND RKGOMfENDATIONS 


The goals of this research program were to: 

• Identify key design parameters that limit high- 
angle-of-attack maneuverability for contemporary 
high-performance attack and fighter-type air¬ 
craft. 

• Postulate fundamental aerodynamic and control 
system design methodologies that will alleviate 
the limiting conditions. 

• Formulate handling qualities requirements for 
high-angle-of-attack maneuvering flight to be 
incorporated in MIL-F-8785C, "Flying Qualities of 
Piloted Airplanes." 

It should be borne in mind that the results are based on variations on a 
single nonlinear aerodynamic model representing a region dominated by 
phenomena that ore highly configuration-dependent. Therefore, the gen¬ 
eralizations drawn are In need of further substantiation. 

Km KEY MANEUVER-LIMITING FACTORS 


The most crucial factors pertain to limiting the safe flight enve¬ 
lope, i.e., departure from controlled flight. By varying static aerody¬ 
namic coupling and cross-coupling coefficients, the nature of high AOA 
characteristics and susceptibility to departure were changed to reflect 
nose slice, rolling, and pitch-up departures. 


Nose slice (yaw) departures predominated for configurations having 
Cngdyn large, negative. Analysis and prior simulation have shown that 
nose slice can be aggravated once sideslip startB by the static aerody¬ 
namic coefficients C£ a and C^. These aerodynamic coupling terms can 

cause a cross-coupled RHP zero to appear in the pitch numerator, n! . 

u stab 

Any simultaneous attempt to control pitch attitude then produces or fur¬ 
ther aggravates directional divergence, 
figurations were prone to flat spins. 


The large negative C na con- 
B dyn 
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Rolling departures predominated when a large component made 

^ n Sdyn positive. These configurations also were prone to wing rock and 
lateral PIO. They exhibited at least two different PSG and spin modes, 
and were sensitive to maneuver and pilot activity at departure onset. 

Pitch due to sideslip, C m g, of either positive or negative sign 
can adversely influence departure and handling characteristics. Posi¬ 
tive C m p can produce pitch-up which tends to mask a "g-break" stall 
warning cue and oppose pilot attempts to prevent stall/departure via 
Torward stick. Negative C m ^ tends “o augment "g-break" stall warning 
but reinforces pilot-commanded recovery pitching moment to the point 
that pitch-roll-yaw coupling may produce violent PSG. The presence of 
appreciable C ra ^ can be detected from flight traces of motion about a 
nominal zero Bideslip via a longitudinal oscillation at twice the dutch 
roll frequency. 

Unstable open-loop characteristic modes (poles) alone were not found 
to limit the safe flight envelope. If the pertinent control numerator 
zeros were favorably located (LHP), the pilot could prevent divergence. 
If the zeros were unfavorably located (RHP) divergence could not be pre¬ 
vented; in this case the presence of open-loop instability increased 
departure severity. Thus, pole-zero locations had influence on depar¬ 
ture susceptibility and severity. In turn, pole-zero locations were 
shifted by static aerodynamic cross-coupling derivatives ' a , N^, Mg) 
when the magnitudes of these coefficients approached those of the common 
static stability derivatives < 4 ;, Ng. Mq): 

• X N^, and Mg, in combination, influence the 
high-frequency modes of both the lateral and 
longitudinal denominators. 
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B. ADDITIONAL FINDINGS AND CONCLUSIONS 


One major contribution has been identification and validation that 
pilot perception of lateral-directional departure susceptibility is 
related to one zero of the numerator ^ becoming negative. Root 
magnitudes more negative than -0.5 rad/sec were consistently rated as 
departure-susceptible, while those less negative (or positive) are rated 
as departure-resistant. This criterion reflects a closed-loop diver¬ 
gence rate limit related to the pilot's threshold for uncommanded motion 
or ability to cope. As such it is a pilot-centered criterion which 
should be applicable for any flight situation, although it has been 
identified in a low-Mach-number, fixed-base simulation. It is consis¬ 
tent with the empirically established airframe-alone departure/spin 
criterion boundaries of Waisaman and extends applicability of that 
criterion to highly augmented airframe cases. It is also consistent 
with previous in-flight simulation of maximum controllable aperiodic 
divergence rates. Finally, it serves as both a design guide and a 
flying quality specification item. 

The "unstable" (RHP) zero generally results from negative (adverse) 
Cng and Cn^g^. It therefore can be alleviated by aircraft configura¬ 
tion and flight control system modifications which principally alter 
those coefficients. As noted above, the zero also can be highly influ¬ 
enced by the static aerodynamic cross-coupling coefficients and Cng. 

A second major contribution has been the development and partial 
validation of a flying quality rating scale applicable to the stall/ 
departure/recovery flight regime. It is believed that the scale has 
potential for defining the aircraft permissible flight envelope in much 
the same manner as the Cooper-Harper scale is now used in defining 
flying quality bounds. 

The simulation demonstrated the Importance of loss-of-control 
warning in pilot assessment of high AOA flying qualities and departure 
susceptibility. A gradual and consistent degradation in stability and 
controllability was found to provide the best warning; however, the 
realization of such warning can be highly dependent upon pitch control 


■i 
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power. Oscillations such as wing rock or lateral PIO cannot be counted 
on for departure warning in air combat if the aircraft can be pulled 
through the '•warning" region before the oscillation becomes apparent to 
the pilot. 

Low-frequency dutch roll or lateral phugoid oscillation was found to 
mask departure onset. Roll about the velocity vector at high AOA pro¬ 
duces significant "nose" yaw. The pilot has considerable difficulty 
differentiating between nose slice onset and low-frequency roll oscilla¬ 
tion about the velocity vector. High roll rate damping (natural or 
augmented) was also found to reduce high AOA departure warning and 
Increase flat spin tendency. 

A flight control or augmentation system mechanization which improves 
and extends the aircraft controllability range by removing natural 
warning can actually be a detriment if the aircraft can still reach some 
departure AOA. To be of real benefit such a flight control system must 
also prevent departure. 

A stick-to-rudder interconnect (SRI) can be beneficial in reducing 
adverse "aileron" yaw, thereby providing more favorable LCDP or 1/T^ 
values. However, if static aerodynamic cross-coupling is strong, the 
analysis of Configuration A 2 has shown that even small sideslip can 
shift the numerator roots in a manner to negate the SRI contribution. 
These results also imply that a finite-authority [3-limiter may be of 
little benefit as a departure preventer if one or more RHP zeros occur 
in control loops the pilot is closing. That Is, pilot control can drive 
the vehicle into a divergence which may exceed the (S-liraiter's capa¬ 
bility. 

Reaults of this study demonstrate that linear or quasi-linear 
frozen-point analysis is applicable to high AOA situations, but consid¬ 
erable caution must be exercised in the Interpretation of results where 
aerodynamics are strong functions of angles of attack and sideslip. Any 
analysis pertaining to such regions must be based on a 6 DOF model with 
aerodynamic and kinematic cross-coupling terms, must include non-zero- 
sideslip trim points, must view transfer-function numerator as well as 
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characteristic-equation parameters, and must be supported by nonlinear 
simulations. 

In particular, the results demonstrated the power of vector polygon 
and closed-loop root extraction analytic techniques for identifying 
cause-effect relationships in highly coupled dynamic systems. The tech¬ 
niques facilitated development of transfer-function literal approximate 
factors for the F-4 in high AO A, non-zero-sideslip flight conditions. 

C. PROPOSED MDL-F-8785C MODIFICATIOHS 

Based upon results of this program it is recommended that modifica¬ 
tions to MIL-F-8785C requirements be considered in three areas. 

1. Departure from Controlled Flight 

The current paragraph (3.4.2.2.1} simply requires that "all Classes 
of aircraft be extremely resistant to departure from controlled flight, 
post-stall gyrations and spina....The aircraft shall exhibit no uncom¬ 
manded motion which cannot be arrested promptly by simple application of 
pilot control." 

The "uncommanded motioi" requirement should be strengthened to 
require that, for roll control input within the service flight envelope 
(SFE) the following stability axis parameter values shall be no more 
negative than: 


Unaugmented airframe: LCDP 





> - 0.01 


Augmented airframe: 1/T^ > -0.5 

An alternate requirement might be that the aircraft should exhibit 
no aperiodic uncommanded motion which exceeds 20 deg/sec and cannot be 
arrested promptly by simple application of pilot control. [The value 
selected here is baaed on a rough average of the simulation pilots' 
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commentary as to their definitions of departure, the spiral divergence 
mode limit of Reference 40, and the 1/T^ limit above.] 

Since the SFE specifically excludes Btall and departure, the above 
criteria would automatically cover high AOA, pre-stall flight. Other 
potential locations for such wording are paragraphs 3.4.2.1.1, Stall 
Approach, and 3.4.2.1.2, Stall Characteristics. 

2. Stall Definition 

The current paragraph (3.4.2.1) allows definition of stall to be 
based on abrupt uncontrollable pitching, rolling, or yawing} or 

intolerable buffet. The results of our piloted simulation indicate that 
any abrupt aperiodic rolling or yawing motion which occurs without being 
preceded by noticeable "g-break" is considered to be a departure, not 
a stall, and results in severely downgraded flying qualities. Thus, 
abrupt roll or yaw motion should be deleted as an allowable definition 
of stall and should not occur prior to stall. If such characteristics 
cannot be achieved with the airframe alone, then the flight control 
system should prevent reaching the AOA at which the abrupt rolling or 
yawing motion is obtained. 

3. Departure Warning 

The current paragraphs (3.4.1.1, 3.4.2.1.1.1, and 3.4.2.1.1.2) 

covering warning pertain mainly to stall. But, since abrupt uncommanded 
roll or yaw is also considered to define stall, the Paras. 3.4.2.1.1.1 
and 3.4.2.1.1.2 requirements, in effect, define departure warning. 
These requirements are based upon fixed margins of V Hta jj or ^Lgtall ant * 
do not take into account the severity of any post-stall/departure mo¬ 
tion, uncontrolled altitude loss, or mission phase. 

Results of this simulation indicated that stall/departure warning 
margin should be related to the severity of any uncommanded motion or 
PSG. Time and altitude loss prior to recovery were heavily dependent 
upon the initial departure severity; pilot commentary indicated the 
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overall departure ratings obtained were heavily influenced by alti¬ 
tude loss and mission phase. Quite obviously, consequences of stall/ 
departure in a ground attack phase can be more catastrophic than in 
high-altitude air combat. Unfortunately, insufficieR'; data ware ob¬ 
tained to establish a specific departure warning crlf ?>'lon. 

4. Compliance Demonstration 

The current specification (Para. 4.1) allows compliance with all 
requirements of Section 3 to be demonstrated through analysis. Although 
compliance with many of the requirements will be demonstrated by simula¬ 
tion, test, or both, this analysis and simulation program has demon¬ 
strated that linearised, frosen-point analytic results may be' midleading 
or extemely difficult to interpret for flight regions where the airframe 
may exhibit large static aerodynamic cross-coupling momenta. Therefore, 
it is recommended that the requirement be changed so that compliance 
with paragraphs pertaining to atall/departure must be demonstrated via 
simulation and/or flight. 
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